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,  pseudorandom '  (PN)^ref erence  code.  Third,  the  informational  signal  may  or 
may  not  be  error  correction  coded  and  hence  requires  an  appropriate 
decoding  process.  , 

*/•.  t  >  -  r*  * 

In  this  study -we  investigateJAthe  utilization  of  a  white-light  or  incoheren 
optical  correlator  for  the  signal  decoding. ^  There  are  several  reasons  for 
selecting  the  optical  technique  rather  than  electronic  or  digital  counter¬ 
part,  because  the  optical  technique  rather  than  electronic  or  digital 
counterpart,  because  the  optical  technique  offers  the  following  capabiliti 
1.  It  is  capable  of  handling  very  long  code  words;  2.  It  has  a  very 
large  space-bandwidth  product;  3.  It  has  the  ability  of  performing 
parallel  processing;  and  4.  The  compatibility  with  optical  storage.  The 
advantages  of  using  white-light  processing  technique  are:  1.  It  elimlnat 
the  coherent  noise;  2.  The  incoherent  system  is  economical  and  easy  to 
maintain;  3.  It  is  suitable  for  color  signal  decoding. 

The  coherence  requirement  of  the  white-light  optical  correlator  for  the 
spread  spectrum  decoding  is  analyzed.  We  have  shown  that  the  temporal 
coherence  requirement  is  strongly  dependent  on  the  spatial  frequency  and 
the  extension  of  the  target  (i.e.,  code  word).  The  spatial  coherence 
requirement  however  depends  only  on  the  extension  of  the  code  word. 

The  effect  on  noise  performance  of  an  incoherent  optical  correlator  is 
also  discussed.  The  output  signal-to-noise  (SNR)  is  evaluated  by  partial 
coherence  theory.  The  sources  of  noise  considered  are  the  grain-noise  and 
the  phase-noise  at  the  input  and  Fourier  planes.  Except  for  the  grain- 
noise  at  the  input  plane,  we  have  shown  that  the  output  SNR  can  be  improve 
considerably  by  using  a  broadband  light  source.  The  overall  SNR  increases 
rapidly  as  the  number  of  spectral  band  filters  increases.  This  increase 
in  SNR  confirms  this  validity  of  noise-suppression  under  a  white-light 
illumination. 

We  have  also  reported  a  polychromatic  correlation  detection  technique. 

This  technique  offers  true  color  correlation  detection  which  is  a  step 
closer  to  actual  human  observation.  The  polychromatic  correlation  may 
also  apply  to  decode  wideband  color  coded  pseudorandom  signal. 

Several  experimental  demonstrations  of  the  incoherent  optical  correlator 
as  applied  to  character  recognition,  pseudorandom  code  are  provided. 
Comparison  with  the  results  obtained  by  the  coherent  system  are  also  given 
We  note  that  the  results  obtained  with  the  incoherent  source  provide  a 
higher  SNR.  Several  demonstrations  for  polychromatic  correlation  detectio 
are  also  included. 

Finally,  the  net  effect  of  this  study  is  to  emphasize  the  truth  of  our 
assumption  that  the  incoherent  or  white-light  optical  correlator  can  be  an 
essential  decoder  for  the  application  of  spread  spectrum  signal  decoding. 
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[.  Introduction 

1. 1  Air  Force  Goals  and  Needs 

The  U.S.  Air  Force  needs  to  improve  the  space  communication  tech¬ 
nique  not  only  by  increasing  the  bandwidth  (e.g.,  spread  spectrum  or 
frequency  hopping)  of  the  transmitted  signal  but  also  needs  to  improve 
the  security  such  that  the  transmitted  signal  is  capable  of  against 
jamming,  fading,  doppler  shifts,  channel  noise,  multipath  effects,  etc. 

The  basic  communication  technique  that  we  would  like  to  investigate 
may  center  in  the  spread  spectrum  communication  technique;  particularly 
in  the  application  of  white-light  optical  signal  processing  to  detect 
a  wideband  pseudo  random  code,  which  was  originally  transmitted  with 
spread  spectrum  communication  method. 

There  are  several  reasons  for  selecting  the  optical  signal  processing 
technique  rather  than  their  electronic  or  digital  counterpart,  because 
the  optical  technique  has  the  following  capabilities: 

1.  It  Is  capable  of  handling  very  long  code  words; 

2.  It  has  a  very  large  space-bandwidth  product; 

3.  It  has  the  capability  of  performing  parallel  information 
processing,  and 

4.  The  compatibility  with  optical  storage  techniques. 

We  would,  however,  use  the  white-light  processing  technique  rather 
than  the  coherent  optical  counterpart,  because  of  the  following  advantages 

1.  It  eliminates  the  coherent  artifact  noise  that  generally 
plagues  the  coherent  optical  system; 

2.  The  white-light  processing  system  is  generally  economical 
to  operate  and  easy  to  maintain. 

3.  The  white-light  technique  is  very  versatile  and  simpler 


to  handle. 
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We  emphasize  again,  the  basic  objective  of  this  research  is  to 
Improve  the  decoding  (i.e.,  detection)  capability  in  such  a  way 
that  high  signal-to-noise  ratio  of  the  correlation  can  be  obtained. 


1.2  Spread  Spectrum  Communication 

Literally,  spread  spectrum  communication  is  one  in  which  the 
transmitted  signal  is  spread  over  a  wide  frequency  band.  Usually 
the  bandwidth  is  much  larger  than  the  base  bandwidth  of  the  informatio- 
signal.  The  most  familiar  example  of  spread  spectrum  modulation  is 
the  frequency  or  phase  modulation.  The  bandwidth  requirement  of  the 
transmitted  signal  is  generally  much  higher  than  bandwidth  of  the 
information  signal.  Following  are  three  general  types  of  spread 
spectrum  modulations^: 

1.  Direct  Sequence  -  modulation  of  a  carried  by  a  digital 
code  sequence  (i.e.,  PN  code)  whose  bit  rate  is  much 
higher  than  the  bandwidth  of  the  informational  signal. 

2.  Frequency  hoppings  -  carrier  frequency  shifting  in 
discrete  manner  followed  by  a  code  sequence. 

3.  Chirp  Modulation  -  the  carrier  frequency  is  swept  over 
a  wide  band  during  a  given  pulse  interval  followed 

by  a  code  sequence. 

Nevertheless,  a  spread  spectrum  communication  system  must  have 
transmitted  bandwidth  much  greater  than  the  information  bandwidth. 

The  art  of  spreading  the  transmitted  bandwidth  is  to  allow  the  com¬ 
munication  system  to  deliver  error-free  information  through  a  noisy 
environment  with  a  proper  decoding  process. 


One  type  of  the  spread  spectrum  communication  that  the  U.S.  Air 
Force  is  particularly  interested  in  is  the  hybrid  forms  of  frequency 
hopping  and  direct  sequence  modulation.  The  form  of  spread  spectrum 
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signal  usually  required  is  a  large  space-bandwidth  PN  code  to  encode 
and  to  decode  the  transmitted  signal.  It  is  this  reason  that  we  believed 
the  optical  processing  technique  would  be  most  applicable  for  a  large 
space-bandwidth  signal.  In  the  optical  signal  processing  technique, 
we  shall  deal  with  an  incoherent  optical  correlator  that  is  capable 
of  performing  the  task  of  correlating  the  PN  code  for  the  spread 
spectrum  signal  decoding. 

1. 3  Report  Outline 

We  deal  primarily  on  an  incoherent  optical  correlator  as  applied 
to  PN  code  correlation  for  spread  spectrum  signal  decoding.  Since 
the  optical  correlation  operation  rely  on  complex  amplitude  operation, 
the  partial  coherence  requirement  of  bhe  correlator  should  be 
evaluated.  The  system  noise  performance  is  also  a  criti'al  issue 
for  the  incoherent  optical  signal  correlator  as  applied  to  the 
proposed  research  problem.  We  shall  report  in  the  fallowing  the  analysis 
of  the  white-light  optical  correlator,  and  the  applications  to  complex 
signal  detection,  pseudorandom  code  (i.e. ,  PN  code)  correlation,  poly¬ 
chromatic  signal  detection  and  noise  performances.  Experimental 
demonstrations  are  also  provided. 


II.  Complex  Signal  Detection  by  Spectrally  Broadband  Source 
2.1  System  Analysis 

We  shall  now  describe  a  technique  of  complex  spatial  filtering  to  be 

carried  out  by  a  spectrally  incoherent  light  source,  for  example,  a  white 
2-4 

light  source  .  We  place  a  diffraction  grating  behind  an  input-signal 
transparency,  s(x,y),  at  the  input  plane  P^  of  a  white  light  optical  pro¬ 
cessor,  as  shown  in  Fig.  1.  The  resultant  complex  amplitude  transmitta.ee 
of  the  input  plane  is  then 

s(x,y)  T(x)  -  K  s(x,y)  [1  +  cos(pox)]. 


(1) 


high  efficient  diffraction  grating 
ex  spatial  filters. 
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where  T(x)  represents  the  diffraction  grating,  K  is  a  proportionality 
constant,  pQ  is  the  angular  spatial  frequency  of  the  diffraction  grating, 
and  (x,y)  is  the  spatial  coordinate  system.  We  note  that,  to  obtain  a 

4 

high  diffraction  efficiency,  a  phase  grating  may  be  used.  Since  the 
input  plane  is  illuminated  by  collimated  white  light,  the  complex  light 
distribution  for  a  given  wavelength  X  at  the  back  focal  plane  P2  of  the 
achromatic  transform  lens  may  be  written 

E(p»q;M  *  CSS  s(x.y)  tl  +  cos(pQx))  exp[-i(px  +  qy)]  dxdy,  (2) 
where  the  integral  is  over  the  spatial  domain  and  spectral  bandwidth  of 
the  light  source,  (p,q)  denotes  the  angular  spatial  frequency  coordinate 
system,  and  C  is  a  complex  constant.  For  simplicity  of  analysis,  we 
dropped  the  proportionality  constant  and  Eq.  (2)  is  therefore 

E(p,q;X)  =  S(p,q)  +  S(p-pQ,q)  +  S(p+pQ,q),  (3) 

mm  mm 

where  s(p,q)  is  the  Fourier  spectrum  of  s(x,y),  p  =  and  q  =  ^0,  (a,0) 
is  the  linear  spatial  coordinate  system  of  (p,q)  and  f  is  the  focal  length 
of  the  achromatic  transform  lens.  In  terms  of  the  spatial  coordinates  a 
and  0,  Eq.  (3)  can  be  written, 

E(cr,P;A)  =  CiSCor.P)  +  C2S(or  -  ^  po,0)  +  C3S(cr  +  ^  P(j,  0).  (4) 

From  the  above  equation,  we  see  that  two  first-order  signal  spectral  bands 
(i.e.,  second  and  third  terms)  are  dispersed  into  rainbow  coloi?  Fourier 
spectrums  along  the  a  axis,  and  each  spectrum  is  centerd  at  a  =  ±(Xf/2n)pQ. 

In  the  analysis,  we  assume  that  a  sequence  of  complex  spatial  filters 
of  various  XQ  are  available ,  i.e., 

H(pn»qn)  =  Kj.  +  K2  +  K  |S(pQ,qQ)|  cos(Poqn  +  <f(pQ,qn)], 

n  =  1,2 . N,  (5) 


where  p  =  fr*’  %  =  X~f^’  ^  t*ie  spatial  frequency  in  the  qQ  direction, 
n1  n1 

K’s  are  the  appropriate  constants,  and  SCp^.q^)  is  the  complex  signal 


spectrum  of  s(x,y).  We  note  that  these  complex  spatial  filters  can  be 

.  6 

computer  generated. 


6 


If  we  place  these  complex  spatial  filters  in  the  spatial  frequency 
plane,  each  centered  at  o  =  (XQf/2n)po,  as  shown  in  Fig.  2,  then  the  com¬ 
plex  light  field  behind  the  spatial-frequency  plane  is 


E(p,q;M  =  S(p-po,q)  2  H(pQ-po ,q^) , 

n-1 


The  corresponding  complex  light  distribution  at  the  output  plane  P3 
of  the  processor  is  therefore, 


g(x,y)  =  2  JT/S(p-po,q)  H^-Po*^)  expfiCp^+q^y) ]dpndqnd\,  (7) 

n=l 


where  the  integration  is  over  the  spatial  domain  and  the  wavelength  spread 
of  the  filters.  We  assume  that  the  signal  spectrum  is  spatial  frequency 
limited  and  the  bandwidth  of  H(p  .q^  is  extended  to  this  limit,  i.e., 

H(pn,q_),  aiSateg 

it  /  _  \  f  nn  s0\ 


H<Pn>V  *  t 


0,  otherwise 


where  =  (X^ f/2n)(pQ  +  Ap),  and  a2  =  (XQf/2n)(pQ  -  Ap)  are  the  upper 
and  the  lower  spatial  limits  of  H(pn,qQ),  and  Ap  is  the  bandwidth  of  the 
input  signal.  The  limiting  wavelengths  of  the  dispersed  spectrums  at  the 
upper  and  the  lower  edges  of  the  filters  are 


p  +Ap  p  -Ap 

X  =  X  ,  and  X,  =  X„  r- 

i  n  pQ-Ap  n  n  pQ+Ap, 


and  the  corresponding  wavelength  spread  over  the  filters  is  therefore, 


AX  =  X 
n  n 


*P0AP 


If  the  spatial  frequency  pQ  of  the  grating  is  high,  then  the  wavelength 
spreads  over  the  filters  can  be  approximated, 


^  ~  — — ^X  ,  p  >>  Ap. 

n  p  n’  ro  ^ 
ro 


(ID 
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Since  the  complex  spatial  filterings  take  place  in  discrete  Fourier  spec¬ 
tral  bands  of  the  light  source,  the  filtered  signals  are  mutually  incoher¬ 
ent,  therefore  the  output  complex  light  distribution  is 

N  ip  x  ip  x  ip  x  ip  x 

g(x,y)~  I  AX  Is  (x,y)e  +  s  (x,y)e  *  s  (x,y)e  0  *  s*(x,y)e  0 

__  i  r  n  u  n  n 

n-l 

ip  x  ip  x  ip  x  ip  x 

+  sn(x,y)e  *  sQ(x,y+Po)e  +  sQ(x,y)e  0  *  s*(-x,-y+PQ)e  °  ], 

(12) 

where  *  denotes  the  convolution  operation  and  the  superscript  *  denotes 
the  complex  conjugate.  For  simplicity,  the  proportionality  constants  K 
were  dropped  from  the  above  equation.  From  Eq.  (12)  we  see  that  the  first 
and  second  terms  represent  the  zero  order  terms,  which  are  diffracted  in 
the  neighborhood  at  (0,0)  in  the  output  plane,  and  the  third  and  fourth 
terms  are  the  convolution  and  correlation  terms,  which  are  diffracted  in 
the  neighborhood  of  (0,-Pq)  and  (0,Pq)  respectively,  as  shown  in  Fig.  3. 
Furthermore,  we  note  that  the  diffracted  output  signal  is  formed  by  inco¬ 
herent  addition  of  the  discrete  spectral  bands,  therefore  the  annoying 
coherent  artifact  noise  can  be  avoided. 

Let  us  now  discuss  the  correlation  term  of  Eq.  (12),  i.e., 

N  . ,  X  f  -iP  y  _ 

R(x,y)  =  I  fJS  S(a  -  2iP0»P)S*(a - fnpo’^e  °  exPU}pf (ax+py) ]dadpdX. 

n=l  n 

(13) 

From  this  equation,  we  see  that  there  is  a  mismatch  in  location  and  in 
scaling  of  the  incoming  signal  spectrum  with  respect  to  the  filter  func¬ 
tion.  In  other  words,  if  the  spatial  carrier  frequency  pQ  of  the  diffrac¬ 
tion  grating  is  high,  a  narrower  spread  in  the  x  direction  of  the  correla¬ 
tion  peak  can  be  obtained.  Thus  the  accuracy  of  the  complex  filtering  in 


the  x  direction  is  somewhat  lower  than  that  in  the  y  direction.  In  other 
words  this  white  light  processing  technique  is  effective  only  in  one  di¬ 
rection,  and  for  some  two-dimensional  processing,  this  technique  may  pose 
some  drawbacks.  However,  this  technique  offers  no  coherent  artifact  noise 
and  it  can  be  applied  to  various  problems  in  optical  information  processing 
We  shall  emphasize  that,  although  the  technique  uses  a  temporal  Incoherent 
source,  the  signal  spectrum  is  displayed  in  partial  coherence  mode  so 
that  the  signal  can  be  processed  in  complex  amplitude. 

We  would  note  that,  the  white-light  or  Incoherent  optical  processing 
scheme  is  to  apply  in  the  detection  of  the  pseudo  random  PN  code  in  the 
spread  spectrum  communication  system.  A  general  matched  filter  as  pro¬ 
posed  by  Caulfield**  may  be  employed  for  this  signal  decoding.  We  shall 
also  note  that,  besides  this  for  the  application  to  decode  the  pseudo 
random  code  for  spread  spectrum  communication,  there  are  several  other 
areas  of  intense  interest  to  the  U.S.  Government  relying  heavily  on  com¬ 
plex  signal  detection  and  hence  stands  to  benefit  from  our  proposed  work. 
Those  areas  are  radar,  sonar,  target  recognition,  etc. 

III.  Coherence  Requirement 
3.1  Basic  Formulation 

Recent  Investigations  of  Morris  and  George^  have  led  to  a  relation¬ 
ship  between  the  correlation  intensity  and  the  wavelength  spread  of  the 

g 

light  source.  Prior  to  their  work,  Watrasiewiez  had  evaluated  the 
effect  of  spatial  coherence  on  the  correlation  intensity  with  a  rectangular 
function.  Here,  we  shall  quantitatively  evaluate  the  temporal  and  spatial 
coherence  requirement  of  a  partially  coherent  correlation  detector,  that 
we  proposed  in  previous  sections. 


s 


The  partially  coherent  optical  signal  processor  of  Fig.  4  employs 

an  extended  broadband  light  source.  For  such  an  optical  processor,  let 

the  source  intensity  distribution,  which  can  be  modulated  by  a  source 
9 

encoding  mask  ,  at  source  plane  PQ  be  y(xo»yo).  In  complex  signal  detec¬ 
tion,  an  input  signal  transparency  of  complex  amplitude  distribution 
t(u,v)  is  inserted  at  the  input  plane  P^  of  the  optical  processor.  A 
matched  spatial  filter  having  an  amplitude  transmittance  H(x,y)  “  T*(x,y) 
is  inserted  at  the  spatial  frequency  plane  P^,  where  T(x,y)  is  the 
corresponding  Fourier  transform  of  the  input  signal  t(u,v)  and  *  denotes 
the  complex  conjugate  from  the  theory  of  partial  coherence^®* the  output 
correlation  intensity  can  be  written  as 

00 

I(u',v')  =  //  Y (x0»yQ)  /  9CA)  c(X)  //  l/x2f2  T(x  +  XQ,  y  +  y^X) 

AS  AA  -® 

H(x,y;A  )exp[i  (xu'  +  yv’)]  dxdy  2  dx  dy  dA,  (14) 

O  AI  O  O 

where  AS  is  the  source  size,  AA  is  the  spectral  bandwidth  of  the  light 

source,  S(A)  is  the  relative  spectral  intensity  of  the  source,  C(A)  is 

the  relative  spectral  response  sensitivity  of  the  detector  or  recording 

material,  Aq  is  the  central  wavelength  of  the  light  source,  and 

H(x,y;A  )  is  the  complex  spatial  filter  at  A  *  A  . 
o  o 

3.2  Temporal  Coherence  Requirement 


We  shall  now  discuss  the  temporal  coherence  requirement  of  the 
light  source.  Let  us  assume  we  have  a  point  source  with  a  broad  spectrum 


i.e.,  Y(xo,yo)  ■  6(xo,yQ).  For  simplicity  we  let  S ( A )  and  C(A)  be  constant. 
The  output  intensity  distribution  of  Eq.  (14)  becomes 


Ku'.v') 


A  +  AA/2 

|A(u’,v’;A)|  dA, 

JX  0  AA/2 


where 


A(u\v';A)  -  If  T(f  ,f  )  T  (f'.f*)  «*p[i2ir(f  u*  +  f  v')]  df  df  . 

x  y  x,  y  x  y  x  y 


*  t2-  t(u’,v')  *  t  Cu*,v*), 

A  O 


(16) 


where  f  ■  x/Af,  f  -  y/Af,  f°  *  x/A  f,  f°  -  y/A  f,  f  is  the  focal  length 
of  the  achromatic  transform  lens  and  *  denotes  the  correlation  operation. 

As  a  special  example  to  illustrate  the  temporal  coherence  require¬ 
ment,  let  us  consider  a  one-dimensional  spatially  Gaussian  target,  i.e., 

2  2 

t(u)  »  e“a  u  ,  (17) 

where  a  is  an  arbitrary  constant.  The  corresponding  Fourier  spectrum  at 
the  spatial  frequency  plane  is 

2, 


T  (f  )  -  /T/a  exp  [-Or  f  )  ], 

X  a  A 

where  f^  *  x/Af,  and  A  is  the  wavelength  of  the  broadband  source  and 
f  is  the  focal  length  of  the  transform  lens.  We  now  insert  a  matched 
filter  for  A  ■  \q  at  the  spatial  frequency  plane. 


(18) 


H(f")  =  T  (f“) 

X  x 


ir/a  exp[-(w/a  f°)  ]. 


(19) 


The  output  intensity  distribution  of  Eq.  (2)  becomes. 


A  +  AA/2  .  2  ,2 

r  0  A  „  A 


I(u’)  *  2ir/a 


2  y  exp  [-2a2  (— ■ -2-  ■ -•)  (u» )  2]dA . 


A  +  A 


(20) 


0  0  0 

From  this  equation,  we  see  that  the  output  intensity  has  again  a 


1(0)  -  2AQ[tan_1(l  +  AA/2Aq)  -  tan_1(l)]. 


Gaussian-like  distribution.  The  correlation  peak  occurs  at  u  *  0,  i.e., 

(21) 

Since  the  spatially  Gaussian  target  essentially  has  low  spatial 
frequency  content,  the  normalized  correlation  peak  (i.e.,  I(0)/AA)  is 
relatively  independent  of  the  spectral  bandwidth,  as  shown  in  Fig.  5. 

To  investigate  the  effect  of  spatial  frequency  on  the  temporal 


coherence  requirement,  we  shall  now  assume  that  a  one-dimensional 


SPECTRAL  BANDWIDTH  AX(A) 

Normalized  correlation  peak  as  a  func 
some  spectral  bandwidth,  for  a  Gauss! 
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sinusoidal  function  as  the  input  object,  i.e. , 


t(u) 


rect(^)  {1  +  co8(2irqu)]. 


where  q  is  the  spatial  frequency  of  the  input  sinusoid,  and 

4  i.  M  <  m. 

rec  Sr  0,  otherwise. 


(22) 


(23) 


where  W  is  the  spatial  extension  of  the  target.  By  substituting  Eq.  (22) 
into  Eq.  (16) ,  we  have 

A(u';A)  -  \J\  j  rect  (^>  [1  +  coa(2irq£)J  rect  (£»■—)  (1  + 


cos[2irqo(£  -  u')]}d£  , 


(24) 


where  Wq  *  WA/Aq,  and  qQ  *  qA/AQ.  Thus  at  u*  *  0,  Eq.  (24)  becomes 


sin(irtfq  ) 

A(W,q,A)  *  A  /A  {W  +  - — -  + 


sinQrWq) 

irq 


sin[irW(q-qo)] 
+  2n(q-qo) 


+ 


sin[irW(q+qo)] 
2it  (q+qQ) 


•}. 


(25) 


The  normalized  correlation  peak  can  be  determined  by  the  following 
equation: 


1(0) 

AA 


,  r\  +  AA/2 

L-  !  0 

AA  J 

A  -  AA/ 2 


(w,q;A)dA 


(26) 


o 

Figure  6  shows  the  variation  of  the  normalized  correlation  peak 
as  a  function  of  the  spectral  bandwidth  of  the  light  source,  for  different 
spatial  frequencies  of  target  (i.e.,  input  object)  signal.  From  this 
figure  we  see  that  the  normalized  correlation  peak  monotonically 
decreases  as  the  spectral  bandwidth  of  light  source  increases.  We  also 
see  that  the  normalized  correlation  peak  drops  rather  rapidly  as  the  spatial 
frequency  of  target  increases. 


Ui  UJ  01- 


SPECTRAL  BANDWIDTH  AX(&) 


On  Che  other  hand.  Figure  7  shows  the  variation  of  the  normalized 
correlation  peak  as  a  function  of  the  spatial  bandwidth  AX,  for  various 
target  extension  W.  Again  we  see  that  the  normalized  correlation  peak 
decreases  monotonically  as  AX  increases.  And  the  normalized  peak  decreases 
very  rapidly  as  the  target  extension  W  increases.  Thus  from  Figs.  6  and  7, 
we  see  that  the  normalized  correlation  peak  strongly  depends  on  the 
spatial  frequency  and  the  extension  of  the  target. 

Using  the  half  power  point  (l.e.,  501)  as  a  criterion  of  the 
normalized  correlation  peak,  we  can  determine  the  spectral  width 
AX  requirement  (i.e.,  temporal  coherence  requirement).  The  results 
are  shown  in  Tables  1  and  2. 

From  these  tables,  we  see  that  for  relatively  small  and  lower 
spatial  frequency  target,  a  relatively  broader  spatial  bandwidth 
of  the  light  source  can  be  used.  In  other  words  the  higher  the  spatial 
frequency  of  the  target  or  the  larger  the  extension  of  the  detecting 
signal  is,  the  higher  the  temporal  coherence  requirement  (i.e.,  narrower 
spectral  bandwidth)  of  the  light  source  is  needed. 

Furthermore,  from  Tables  1  and  2,  a  plot  of  the  space  bandwidth  product 
(i.e.,  Wq)  as  a  function  of  the  required  spectral  bandwidth  of  the  light 
source  is  shown  in  Fig.  8.  From  this  figure,  we  see  that  the  space 
bandwidth  product  exponentially  decreases  as  the  spatial  bandwidth  of 
the  light  source  increases.  Thus,  for  a  large  space  bandwidth  product, 
an  extremely  high  temporal  coherent  light  source  is  required. 

3.3  Spatial  Coherence  Requirement 

We  shall  now  discuss  the  spatial  coherence  requirement  of  the  light 
source  for  the  correlation  detection.  For  simplicity,  we  assume  that  the 
light  source  Y(xQ,yo)  is  one-dimensional  extended  monochromatic  source. 
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Figure  7.  Normalized  correlation  peak  as  a  function  source 

spectral  bandwidth,  for  different  target  extensions. 
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y(xo,yo)  =  rect  (— )  f(y<)) 


where 

x  .  1,  lx  I  <  AS/2, 

rect  C^g)  -  {  ( 

0,  ohterwise, 

and  A  is  the  size  of  the  extended  source.  Similarly,  the  output  plane 


intensity  distribution  can  be  shown  as 
AS/ 2 


I(u') 


rectOrf-)  A2(uf;x  )dx  -  2 
AS  Oo 


-AS/ 2 


A2(u';x  )dx  ,  (29) 

o  o 


where 


A(u’;xo) 


t  U)  t(C-u’)exp[i2irfx  £]<*£, 

o 


and  f 


x  X  f  * 
o  o 


We  shall  again  consider  a  spatially  Gaussian  target  of  Eq.  (17), 


where  we  obtain 


2  2  2 

A(u';x  )  ■-  j  exp [-( — - —  f  )  ]  exp [-  j  (u1)  ] 

“  a  A  /  A 

a  /  2  o 


The  corresponding  output  intensity  distribution  can  be  shown  as 

AS/2 


I(u’) 


\exp[-a2  (u’)2] 


exp[-(—  f  )  ]dx  . 
v  a  x  o 

o 


Since  the  integral  is  Independent  of  u',  it  can  be  shown  that  the  output 

intensity,  like  before,  has  a  Gaussian-like  distribution.  To  determine 

the  correlation  peak,  we  shall  let  u'  **  0,  i.e., 

AS/2  2 

1(0)  -  exp [ - (c~  f  )  ]dx  .  (3! 

4  a  x  o 

a  ^ 

It  is  evident  from  the  equation  that  1(0)  depends  on  the  size  of 
the  extended  monochromatic  source  AS.  Figure  9  shows  the  normalized 
correlation  peak  1(0) /AS  of  the  Gaussian  target.  We  see  that  the  cor¬ 


relation  peak  drops  quickly  with  the  increasing  size  of  the  source. 


We  shall  now  consider  the  effect  of  spatial  frequency  in  the  corre¬ 
lation  peak.  Again  we  use  a  spatially  limited  sinusoidal  grating  of 
Eq.  (22)  as  the  input  signal.  For  simplicity  of  notation  we  let 
q  ■  xo/(XQf)  and  g  ■  ir/(AQf),  Eq.  (33)  becomes 


1(0)  -  2 


AS/2 

► 

. 

0 


(0;x3dx  , 
o  o 


(34) 


where 

sin(gWx  )  sin[W(gx  -irq)]  sin[W(gx  +irq)] 

A(0,x  )  -  - —  +  - 2_ - +  - — 2_ - 

o  gx  (gx  -  irq)  (gx  +  irq) 


1  2  qsintW(gxo-2itq)] 

(gx0+2lrq)  «*,)  +  - (gx  -  2 


Trqsin(Wgx  ) 
+  - — ] 


(35) 


q  is  the  spatial  frequency  and  W  is  the  spatial  extension  of  the  target. 


Since  q  is  not  related  to  xq,  the  spatial  coherence  requirement  is 


evidently  independent  of  the  spatial  frequency  of  the  input  object,  as 
expected.  However  the  spatial  coherence  strongly  depends  on  the  extension 
of  the  target,  as  can  be  seen  in  Fig.  10.  From  this  figure,  we  see  that, 
the  normalized  correlation  peak  monotonically  decreases  with  increasing 
source  size.  The  rate  of  the  decrease  rapidly  increases  for  larger 
extension  of  the  target.  Again  by  using  the  half  power  point  criterion, 
the  spatial  coherence  requirement  for  various  values  of  the  target 
extension  are  tabulated  in  Table  3.  The  table  shows  the  spatial  coherence 
requirement  (i.e.,  source  size  requirement)  increases  rapidly  as  the  target 
extension  increases. 

In  concluding  this  section,  we  have  developed  equations  for  the 


coherence  requirements  of  correlation  detection  for  partially  coherent 
optical  processing.  We  have  also  calculated  numerical  estimates  for 


Figure  10.  Normalized  correlation  peak  as  a  function  of 
source  size,  for  different  target  extensions. 
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these  requirements  using  half  power  criterion  for  some  special  cases. 
We  have  shown  that  the  temporal  coherence  requirement  strongly 
depends  on  the  spatial  frequency  and  the  extension  (i.e.,  space 
bandwidth  product)  of  the  target  signal  to  be  detected.  However, 
the  spatial  coherence  requirement  depends  only  upon  the  extension 
of  the  target.  We  have  also  shown  the  coherence  requirement  for 
a  spatially  Gaussian  signal,  as  a  special  case.  For  this  case,  the 
normalized  correlation  peak  is  relatively  independent  of  the  spectral 
bandwidth  of  the  light  source,  but  it  decreases  rapidly  as  the 
source  size  increases. 


IV.  Noise  Performance 
4.1  Problem  Formulation 

The  noise  performance  of  an  achromatic  coherent  optical  system 

12 

was  reported  by  Leith  and  Roth  .  They  approached  the  analysis  through 
the  introduction  of  the  concept  of  a  three^-dimensional  transfer  function 
to  describe  the  noise-suppression  properties  of  the  system.  They  showed 
that  such  a  system  produces  considerable  noise-suppression  if  the  system 
is  illuminated  by  a  broad-spectrum  light  source.  Both  signal -dependent 
and  signal-independent  noise  were  considered.  In  each  case,  they  showed 
that  the  achromatic  coherent  system  behaves  much  like  an  incoherent 
imaging  system. 

Here,  we  shall  analyze  the  noise  performance  for  a  white-light  optical 
1-3 

signal  processor  (or  a  partially  coherent  optical  processor),  as  proposed 
in  Fig.  11.  We  shall  derive  the  output  signal-to-noise  as  a  means  of 
measurement  for  the  noise  performance  of  the  system.  The  noise  sources  that 
we  shall  discuss  are  primarily  due  to  grain-noise  and  plane-noise  at  the 
input  and  Fourier  planes.  The  development  of  the  output  signal  irradiance 
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and  noise  fluctuation  are  primarily  based  upon  the  partial  coherence 
theory  of  Wolf^’^.  Since  output  noise  fluctuation  is  dependent  on 
the  degree  of  coherence  illumination,  we  shall  quantitatively  analyze 
the  effect  of  noise-suppression  due  to  temporally  and  spatially  partial 
coherent  illumination.  In  this  section,  we  shall,  however,  analyze 
the  noise  performance  due  only  to  the  temporally  coherent  illumination. 

The  noise  effect  due  to  spatially  coherent  illumination  will  be 
followed  by  a  subsequent  report. 

We  shall  now  evaluate  the  noise  performance  of  the  proposed  white- 
light  optical  signal  processor,  as  described  in  Fig.  11.  We  shall  first 
investigate  the  statistical  nature  of  intensity  distribution  at  the 
output  plane.  In  order  to  do  so,  we  shall  evaluate  the  output  light  intensity 
distribution  due  to  the  nth  narrow  spectral  band  spatial  filter,  such  aS13 
X^n  oo  os 

1  n  (x .  tf')  « J  JJ  r  i. .  X)SC^)C(x)  |  { $S(x0+*-\fV0 ,  yft) 

Xin.  *  °° 

Hn  (<*,  (0  •* pH  tH***  a°ia(i  I 

(36) 

where  Y(xo>yQ)  is  intensity  distribution  of  the  white-light  source, 

S(\)  and  C(A)  are  the  relative  spectral  intensity  of  the  light  source 
and  the  relative  spectral  response  sensitivity  of  the  detector,  respec¬ 
tively;  A^  and  A  are  the  longest  and  the  shortest  limiting  wavelenths 
of  the  nth  spectral  band  filter,  S(xq  +  a  -  Afv  ,  yQ  +  8)  represents 
the  Fourier  spectrum  of  the  input  signal  s(x,y)  due  to  wavelength 
A,  vQ  in  the  spatial  frequency  of  the  phase  grating,  (a, 3)  is  the 
nth  spectral  band  spatial  filter  and  f  is  the  focal  length  of  the 
achromatic  Fourier  transform  lens. 
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The  overall  intensity  distribution  for  the  entire  spatial  band  of 
the  light  source  would  be, 

IK*')  “I  InK*') 

*t*l  _ 

M  *-An  00  00 

~l  S  tf.  i  V(1) 

**•  \ln .« 

(37) 

where  N  is  the  total  number  of  the  narrow  spectral  band  filters. 

The  overall  output  signal-to-noise  ratio  (SNR)  of  the  white-light 
processor  can  be  written  as 


SNR 


4-1  snr* 


where  SNR  represents  the  output  signal-to-noise  ratio  due  to  the  nth 
n 

spectral  band  spatial  filter. 

We  shall  now  consider  a  spectrally  broadband  point  source  (i.e.,  a 
temporal  incoherent  source)  for  illumination.  The  spatial  intensity 
distribution  of  the  light  source  can  be  described  by  a  5-function 
[i.e.,  Y(xQ,yo)  *  6(xo>yQ)].  For  simplicity,  we  assume  that  C(X) 
and  S(A)  are  constant  over  the  spectral  width  of  the  light  source. 

Thus  Eq.(37)  can  be  written  as 

(39) 

where  the  subscript  n  represents  the  nth  output  irradiance  due  to  nth 
spectral  filter  Hn(a,8)  and  the  proportionality  constant  is  disregarded 
for  convenience.  In  analysis,  we  let  the  input  signal  be  a  one  dimen¬ 
sional  object  of  y  variable  [i.e.,  s(y)].  The  phase  grating  at  the  input 
plane  is  perpendicular  to  y  direction,  which  can  be  described  as 


exp(i2irvox) .  The  corresponding  Fourier  smeared  spectra  are  also  one- 
variable  function  in  the  8-direction  and  they  are  smeared  into  rainbow 
color  in  the  a-direction.  Let  us  denote  that  the  width  of  the' nth 
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narrow  spectral  band  filter  be  Ax^.  If  this  filter  is  placed  in  the 
appropriate  spectral  band  of  the  smeared  Fourier  spectra,  the  corres¬ 
ponding  spectral  bandwidth  AXn  is  obviously  determined  by  the  spatial 
width  Aan  of  the  filter,  which  can  be  approximated  by  the  following 
equation; 

fcXn  =  (X*n-  XJL*)  *=  /v.  5  (40) 

liy  fixing  the  spatiai  width  Act^  oL  the  filter,  the  redundancy  of  the 
filtered  signal  would  be  proportional  to  the  spectral  bandwidth  AX^  of 
the  filter.  In  other  words,  the  broader  the  spectral  bandwidth  AX^, 
the  higher  the  redundancy  of  the  filtered  signal  is  expected.  From 
Eq.(AO),  it  is  evident  that  either  a  lower  spatial  frequency  vq  of 
phase  grating  or  shorter  focal  length  f  of  the  transform  lens  should 
be  used  to  improve  the  output  SNR  at  the  output  plane.  However  de¬ 
creasing  the  grating  spatial  frequency  vq  also  reduces  the  number 
of  spectral  band  filters  in  the  Fourier  plane.  Thus  it  would  also 
reduce  the  processing  capability.  Although  the  output  SNR  would 
decrease  as  increases,  however  larger  the  vQ  is  used  also  implies 
that  a  larger  number  of  spectral  band  filters  can  be  utilized  in  the 
Fourier  plane.  Thus  the  overall  output  SNR  would  be  increased  due 
to  the  addition  of  the  mutually  incoherent  filtered  signals,  from 
each  spectral  band  filters  Hn(8). 

It  may  be  known  that  the  output  intensity  fluctuation  is  primarily 

due  to  the  statistical  nature  of  the  phase  and  the  amplitude  noise 

fluctuations.  We  shall  now  define  the  output  SNR  as  an  ensemble  average 

of  the  image  intensity  I  (y)  divided  by  the  square  root  of  its  variance 

ar  ,  that  is 
In 

SNRn(^)=E(I»Cj'))/<ru 


f 


(41) 


(42) 


^Xn  — j" 

•  • 

Again  we  note  that  a  one  dimensional  notation  is  used,  since  the  input 

signal  is  assumed  to  be  one-dimensional.  In  order  to  describe  the  output 

SNR  explicitly  due  to  the  effect  of  partially  coherent  illumination,  we 

shall  define  a  normalized  signal-to-noise  ratio  SNR,  such  as, 


sN«„(ao  = 


j-OY  pcuir+uJUt^  CpKiJYtnX  wn. 


which  is  equal  to  the  SNR  of  the  partially  coherent  case  divided  by 
the  coherent  case,  where  the  subscript  n  denotes  the  nth  spectral 
band  filtering. 

We  stress  that  the  spatial  average  of  Eq.  (43)  shall  be  used  for 
estimating  the  SNR  over  the  entire  output  image  plane,  for  example; 
(a)  For  a  continuous  representation,  the  output  normalized 
SNR  is , 


v* _ 


-Wz 


> 


(44 


where  L  is  the  size  (i.e.,  the  length)  of  the  output  image,  and 

(b)  For  a  discrete  representation, 

_  . 

SN«n  =  fJ  l  SNR„(tf:)  ,  <45 

m  *=■  | 

where  M  represents  the  total  number  of  sampling  points.  From  these 
definitions  it  is  reasonable  to  assume  that  the  SNR^  is  equal  to  unity 
for  a  strictly  coherent  illumination  and  SNRn  increases  as  the  degree 
of  coherence  decreases. 

Since  SNR^  is  defined  for  nth  spectral  band  filter,  the  overall 


output  normalized  signal-to-noise  ratio  can  be  written  as 
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SNR  =-frZ  sn*«  , 


where  N  is  the  total  number  the  spectral  band  filters  in  the  Fourier 
plane,  which  can  be  estimated  by. 


AX  and  X^  are  the  spectral  bandwidth  and  the  mean  wavelength  of  the  light 

source  (e.g. ,  white-light  source),  respectively,  v  is  the  spatial 

frequency  of  the  phase  grating,  andAv^  is  the  spatial  frequency  bandwidth 

of  the  input  object,  in  a  direction. 

We  stress  that,  Eqs.  (39)  to  (46)  shall  be  used  for  the  evaluation 

of  the  noise  performance  of  a  white-light  optical  signal  processor. 

Prior  going  into  the  evaluation  of  the  SNR,  we  shall  specify 

the  statistical  properties  of  the  phase  noise  and  film-grain  noise, 

14 

respectively.  Let  us  now  adopt  an  overlapping  grain  mode  to., 
describe  the  granularity  of  the  recording  medium.  For  simplicity  of 
evaluation,  we  would  make  the  following  assumptions: 

(a)  The  film-grains  are  rectangular  in  shape  and  the  center 

of  these  rectangle  grains  fall  randomly  over  the  recording  plate,  as 
shown  in  Fig.  12. 

(b)  The  distribution  of  the  grains  (i.u.,  the  choice  of  where 
the  center  of  a  rectangle  falls)  is  assumed  to  be  independently 
distributed. 

and  (c)  The  distribution  of  the  grain,  over  the  recording  plate, 
is  Poisson  distributed,  i.e.. 


Q  (**)  = 


m  !  (  N  -  m)  | 


pVp>"-" 


The  mean  value  and  the  second-order  moment  (correlation  function)  are: 


E  (<*(•(>] 


(4») 


'  Ji  rfp  ■ 

_  4  % 

ja  ml 

fi£3  mT^  ■■ 


*a  f  '  d 

3p-  '  #  V' 

^4  \\v  ! 


I! 
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and 

£[$<*.)<,<*.)]  =C^  e«p[?.3DR(|^,-^|/A)J  ;  (50) 

respectively,  where  D  is  the  photographic  density,  i  is  correlation 
distance  and 


V. 


0 


» 


for  Ui-fxh-8- 

(51) 


where  the  function  R C | y ^  -  y ^ |  /-8-)  is  shown  in  Fig. 13. 

Although  the  film  granularity  is  generally  a  nonstationary  Markow 
Chain  for  simplicity,  we  assume  that  the  grain-noise  is  a 

stationary  process. 

As  for  the  phase-noise,  the  phase  fluctuation  <Ky)  is  primarily 
due  to  the  thickness  variation  and  the  refractive  index  fluctuations 
of  the  emulsion.  For  simplicity,  we  would  make  the  following  assump¬ 
tions: 

(a)  >(y)  is  a  stationary  random  function  of  the  coordinates, 

the  probability  distribution  of  4> (y)  is  identical  throughout  the  entire 
recording  medium. 

(b)  It  is  reasonable  to  assume  that  <Ky)  has  zero  mean  i.e., 

E [ (y )  ]  =*  0,  since  the  phase  fluctuation  can  be  regarded  as  a  bipolar 
function. 

(c)  The  autocorrelation  function  of  <p (y )  is  dependent  unly  on 
the  difference  in  distance  )y^  -  y  |  since  (y)  is  assumed  to  be  a 
stationary  process.  We  note  that  the  autocorrelation  function  has 
various  functional  forms;  for  example  [13j,  the  Gaussian  form  represen¬ 
tation; 


> 


(52) 


ariat 
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and  the  exponential  form  representation; 

E  =o2exp[-|u,-ift|/d.]  ;  <53) 

2 

where  a  is  the  variance  of  $(y),  and  the  quantity  d  is  called  the 
scale  size  of  the  irregularities  and  can  be  regarded  as  the  size  of 
an  average  blob,  in  the  emulsion.  However  d  is  best  known  as  the 
correlation  distance.  For  simplicity,  we  shall  however  use  the  exponen¬ 
tial  form  for  evaluating  the  output  SNR. 

4.2  Effect  of  Noise  at  Input  Plane 

We  shall  now  evaluate  the  noise  performance  of  the  proposed  white- 
light  optical  signal  processor  due  to  input  plane.  Let  us  consider  the 
output  intensity  distribution  due  to  nth  spectral  band  filter  of 
Eq.  (37).  By  Perseval's  Theorem,  the  output  intensity 
distribution,  in  a  one-dimensional  form,  can  be  written  as 

I,CJ7  =  I  I  <54) 

where  the  subscript  n  denotes  the  nth  spectral  band  filter,  *  denotes 

the  convolution  operation  and  t(y’)  represents  the  one-dimensional  input 
signal  with  complex  noise.  The  complex  noise  is  due  to  the  random 
thickness  fluctuation  of  phase  distribution  4>(y)  and  the  granularity 
G(y)  of  the  input  transparency.  We  note  that  the  phase-noise  and  the 
grain-noise  can  be  regarded  as  multiplicative  and  additive  type  noise, 
respectively.  The  overall  input  signal  can  be  written  as, 

•Klj)”  S(fl)exp[i.K<t>(jp]  +  Q(^)  ,  (55) 

where  k  =  2tt/.\,  and  s(y)  is  a  one-dimensional  signal. 

To  evaluate  the  output  SNR,  we  let  the  input  signal  s(y)  be  a 
sinusoidal  grating  with  a  spatial  frequency  v.  Thus  Eq.  (55)  becomes, 
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t(^)  —  (l  +  S^2nj7^)exp(cK4>(y)]  +  CT(p 


(56) 


We  shall  now  Investigate  the  effect  due  to  phase-noise  and  film- 
grain  noise  independently  in  the  following: 

4.2.1  Phase-Noise 

To  investigate  the  effect  due  to  phase-noise  at  the  input  plane,  we 
simply  let  the  film-grain  noise  be  zero  and  assume  that  the  nth  spectral  band 
filter  Hn(B)  be  bandlimited,  i.e., 

H,(<°  =  rect  (-J^T )  . 

By  substituting  Eq.  (57)  into  Eq.  (54)  we  have 

\fclt  oo 

In(!j')=f  |  J  (l  +  St,2«vg)e*p[i.K<f>(J) 

OO 

If  we  take  the  ensemble  average,  Eq.  (58)  becomes 


Mn  oo 

E  (Kf)J  =1  f 

Aju.  -0'7 


(59) 


We  shall  evaluate  Eq.  (56),  with  the  following  assumptions: 

(a)  k  $(y)  <  1,  a  weak  phase-noise,  and 

(b)  higher-order  moments  of  <fr(y)  are  much  smaller  than  the 
second-order  moment. 

With  reference  to  the  Taylor  series,  the  ensemble  of  the  exponen¬ 
tial  term  can  be  written  as: 


E  [exp  [  6  k  $<^313 

=.|  +  CKE [♦<*>- 4>c^a>>  t(c k>* E  ( MW- 4>t*i)f ]  ■  •  •  , 

(60) 


By  utilizing  the  exponential  form  representation  of  Eq.  and 

the  assumption  of  E ( d> (y ) 3  =  0,  Eq.  (57)  is  reduced  to 

E  [exp  [Ck£$(^)-4>(^  )]}]*»  |  -  Kzcrz+  K*<y2exp(-l^-y  |/d)  (6i) 

By  substituting  Eq.  (61)  into  Eq.  (59),  we  have 


E  tup] - i:  ( i)  r o-  ♦  w)i»<pj 


XjUv 
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=axxf  C^'j- 


4tit<r’ 


C\gn+  AA.) 


(62) 


where 


I  Cfl')  Aj(hSi~MVy) 


(62a) 


i,(^”.i(i+suw*)5T^^top(^^<i? , <62t> 


Xnd. 


[,()(')  =  j  (l*Sw2TH>j)  $<~2'U|(fl-9) 


e«P(-K/A)4'j  (62c) 


The  ensemble  of  the  square  irrad lance  can  be  obtained  by  the  following 
integral  equation: 


oo  t 

E [ill <*>]  -  E(  { I  SJ  ( 1 4 aUwvj )(n )S-^gy 
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'J 1  e«p{1  K  [»(?)-  *( f)J}  Aj  A  a  4x)*  3 
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,J$  iX,iXt  Sffi  (nSUJ»)l'J,)(l+Sl»ZirVty)(l*SiU2tlU*|t) 

X  S^CTy ' C*  ~ ?»>  E[<up[iK,t*M.)-»(?,fl 

(63) 


Similarly,  the  following  equations  can  also  be  obtained  for  either 
the  weak  phase-noise  or  the  small  higher  order  moments  assumptions, 
such  as 

E  [<wp(  £  K,  C4>t^.)- +  i  KS  +  CJOJ 1  ] 

=  l-lKNKj)s,-K?E[4(30<K^)]-Kh[4'(i,)<(>(j1)] 

♦  ic,  kj  s  C<l>(a, )<(>(?,)) -e  C*rj.)4>(^)j 


(64) 
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By  substituting  Eqs.  (64)  and  (50)  into  Eq.  (63)  and  using  the  symmetrical 


properties  of  the  integral  function  about  the  variables  (y^.y£)  and 
(y ■  ,  y .  we  obtain  the  following  result: 

2  2  l  * 

E  =JJ  I*(y)-  <y>-Ii(y'>U(y')3 

XjUl 

*  (  -7? - JT ) 

A, 

2 

=  (^»-^i«)*I|'4l^)-  if  <-i)-  <65) 

^An*AjUv 

where  1^  to  1^  are  defined  in  Eq.  (62).  With  reference  to  the  definition 
of  SNRn  of  Eq.  (41) , the  signal-to-noise-ratio  due  to  the  nth  spectral 

E  [l,(m3 _ 

0«  -xzrHV)- 

(66) 


band  filter  would  be 

CO 

SNRn(^')* 


where  the  superscript  (1)  represents  the  output  SNR  due  to  the  phase- 

noise  at  the  input  plane  and  AX  =  A.,  -  A.  is  the  spectral  bandwidth 

n  hn  itn 

of  Hn(B).  With  reference  to  the  definition  of  Eq.  (45),  the  output 
normalized  SNR  can  be  obtained,  such  as 


- (0 

SNRn 


(67) 
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We  note  chat  the  normalized  output  signal-to-noise  ratio  SNR  ^ 

due  to  nth  spectral  band  filter  shall  be  used  to  evaluate  the  input 

phase-noise  performance  as  a  function  of  AA^  (i.e.,  the  temporally 

partial  coherence  effect)  of  the  proposed  white-light  optical 

signal  processor  of  Fig. 11.  We  shall  first  investigate  the  behavior 

of  SNR  ^  as  a  function  of  spectral  width  AX  ,  for  various  values 
n  n 

of  normalized  spatial  frequency  ft,  i.e.. 


where  v  is  the  spatial  frequency  of  the  input  object  and  is  the 

cut-off  spatial  frequency  of  the  optical  processor.  Plots  of  the 

SNR  ^  as  a  function  of  AX  ,  for  various  values  of  ft  are  given  in 
n  n 

Fig. 14.  From  this  figure,  we  notice  that,  the  higher  the  spatial 
frequency  of  the  input  object  the  higher  the  SNR^^1^  for  broader 
bandwidth  AAr  (i.e.,  broader  band  temporally  partial  coherent  illumina¬ 
tion).  In  other  words,  the  SNR  ^  would  not  noticably  improve  for 

n 

low  spatial  frequency  object,  for  phase-noise  effect  at  the  input 

plane.  Nevertheless,  the  situation  would  be  quite  different,  as  we  will 

show  in  a  moment,  for  multi-spectral  band  filterings  as  depicted  in  Fig. 11. 

Let  us  now  investigate  the  SNR  ^  as  a  function  of  AX  due  to 

n  n 

A  (i.e.,  shortest  wavelength  limit  of  the  nth  spectral  band  filter). 

FigurelS  shows  a  set  of  plots  of  SNR  as  a  function  of  the  spectral 

bandwidth  AA  ,  for  various  values  of  X  .  From  this  figure,  we  see 
n  x*  n 

that,  the  SNR  ^  increases  as  the  spectral  bandwidth  AX  increases, 
n  n 

For  a  shorter  wavelength  limit  A.  ,  the  SNR  shows  a  noticable 

*n  n 

improvement  for  wider  spectral  bandwidth  AX^  of  this  filter.  In 
other  words,  the  overall  phase-noise  improvement  is  more  noticeable 


for  short  wavelengths. 
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We  shall  now  investigate  the  effect  of  overall  phase-noise  performance 
of  the  proposed  white-light  optical  processor  utilizing  multispectral  beam 
filterings.  With  reference  to  Fig.  .11,  the  overall  output  signal-to-noise 
ratio  of  the  processor  would  be  the  sum  of  the  mutually  incoherent  narrow 
spectral  band  SNRn.  Thus  the  overall  output  normalized  signal-to-noise 
ratio,  due  to  phase-noise  at  input  plane,  can  be  written  as 


5NRCI>  =  -frI  snr"’ 

w=| 


_ J _ L  y  y  r  -  +  I|  i 


A  plot  of  the  overall  SNR'  as  a  function  of  the  number  of  spectral 
band  filters  N  is  given  in  Fig. 16.  For  this  figure,  we  see  that,  the 
overall  SNR^  increases  rapidly  (approximately  exponential  increasing) 
as  the  number  of  spectral  band  filters  increases.  This  effect  is  quite 
noticeable  for  broad  spectral  band  filtering  as  in  contrast  with  the 


narrow  spectral  band  filters 


Thus  we  conclude  that,  the  broad¬ 


band  white-light  optical  signal  processor  does  improve  the  phase-noise 
performance  at  the  output  plane,  as  in  contrast  with  a  coherent  processor. 
4.2.2  Grain-Noise 

Let  us  now  evaluate  the  effect  of  the  noise  performance  due  to  film- 
grain  noise  at  the  input  plane.  We  note  that,  by  varying  the  wavelength 
of  the  light  source,  it  is  generally  noticeable  to  make  a  distinction 


between  the  grain-noise  and  the  object,  since  the  photographic  grains 
will  not  alter  by  different  wavelength  illumination.  In  other  words, 
a  temporally  partial  coherence  illumination  can  not  reduce  the  effect  of 


film-grain  noise  at  the  input  plane.  Nevertheless,  we  shall  quantitatively 
evaluate  this  effect  in  the  following: 


As  similar  to  the  phase-noise  evaluation,  the  output  intensity  dis 


tribution  can  be  obtained  by  the  following  equation: 


is  the  input  signal  function.  By  taking  the  ensemble  average  of  the 
above  equation,  we  have 

00 

E  I  (  jjjl  +  S2nW^)(l4SW*V^  +  (lt5«<n.*lTVJ) 

■"■(3-3)  ' 


With  reference  to  the  definition  of  Eq.  (50),  we  see  that 
E[£(y)  G(y)]  is  discontinuous  at  (y-y|  >  i.  The  integrant  of  Eq.  (70) 
is  not  easy  to  evaluate. 

To  simplify  the  operation,  we  approximate  E[G(y)  G(y)]  by  a 
continuous  function,  as  shown  below, 

E  feC‘j)C1(^)]«C^{(e*p(2.3D)-l]expC-|^5|A)^]  f  ( 

where  *  E[G(y)],  and  a  =  l/{2.3  +  ln[ exp (2 . 3D) -1 ] }  is  a  continuous 
function. 

To  justify  the  approximation  of  Eq.  (71),  a  plot  of  Eq.  (71)  as 
compared  with  variation  of  Eq.  (50)  is  shown  in  Fig.  17.  From  this 
figure,  we  see  that  only  small  amounts  of  deviation  is  introduced 
by  Eq.  (71)  as  compared  with  the  variation  of  Eq .  (50).  Thus  the 


fc  A 
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approximated  function  of  Eq.  (71)  would  not  expect  to  introduce  significant 
errors  in  the  subsequent  evaluation. 

Now  by  substituting  Eq.  (71)  into  Eq.  (70),  and  using  the  same  approach 
as  Eq.(63),  the  ensemble  average  of  the  output  irradiance  ln(y'),  due 
to  the  nth  spectral  band  filter,  can  be  shown  as 


E  U  n  <  O)  =  if  ( »>  2  ca  l»  ♦  ej  [«p  (md)-  i  ]  V  i,  <  10 


(72) 


where 


.L  "(y-V  "  > 

ViP  =1  . 


(72a) 


(72b) 


and  I^(y ' )  is  given  in  Eq .  (62a).  Similarly,  the  formula  for  the 
covariance  can  be  derived  in  the  following: 

Xln  -oo 

*  (l  +  Si^.2rfV^2)(|  +  S^*2«^^)+A(l45^n^^) 

,Aii  )  < 
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where  we  have  utilized  the  symmetry  property  of  the  integral  function 
with  respect  to  the  variables  (y^j^)  and  (y2>y2).  We  also  assume  that 
the  higher-order  moments  of  the  random  function  C(y)  are  much  smaller 
as  compared  with  the  second-order  moment  (autocorrelation)  function. 
Furthermore,  by  substituting  Eq.  (?1)  into  Eq.  (73),  we  obtain 

E  1 1*  ( *')]  =  if  <  D')  -  4  CA  if  (  f)  I.C j- )  *  6  4  [exp(2.3D)- 1 J 

where  =  E[G(y)]. 

With  reference  to  Eqs.  (41)  to  (45),  the  output  normalized  SNR,  due  to 
the  input  grain-noise  can  be  written  as 

SNRn  “-“T-Z  IX  (</"») 

•V|  w«l 

4  (exp  (2.3 D)  -  I  )l,(^)lx(^)]} K  t4  30) 

-I  ]  I, 

+  [exp(2.3o)-  I ]  J<( J-)Ir(£)]-Cj  [I02( tfi)  * 

(6xp(2.30)-l)l4(^L)lr(^)32j,/2  ^  (75) 

where  the  superscript  (2)  denotes  the  SNR  due  to  the  input  grain-noise. 

We  note  that  if  <<  1,  then  Eq.  (75)  reduces  to 

^*> _ Lr  _ _ 

5  n  Hi,  2Ca  Ue»p  (MUD.  01,  Whitin'*  t  <76) 

From  the  above  equation,  it  is  evident  that  the  SNR  ^  is  independent 
of  the  spectral  bandwidth  AA^  (i.e.,  the  degree  of  temporal  coherence). 


In  other  words,  the  effect  of  the  grain-noise  performance  (due  to  the 

input  grain-noise)  can  not  be  improved  by  using  a  temporally  partial 

coherent  or  simply  temporally  incoherent  illumination,  as  expected. 

We  note  that  this  behavior  is  quite  different  from  that  as  encountered 

with  the  input  phase-noise  case,  as  described  previously.  There  are, 

however,  several  interesting  behaviors  that  can  be  infered  from  the 

-  (2) 

result  of  Eq.  (76).  For  example,  the  SNRn  is  dependent  on  the  mean 

-  (2) 

value  CA  of  the  input  grain-noise.  The  value  of  SNR^  decreases 

(approximately  exponentially)  quite  rapidly  as  the  mean  value  CA 

-  (2) 

increases,  as  shown  in  Fig.  18.  However  the  SNR^  is  independent 

of  the  spectral  bandwidth  of  the  processing  filter,  as  shown  in 

Fig.  19.  From  this  figure,  we  see  that,  the  overall  signal-to-noise 

- (2) 

(i.e.,  SNR  )  should  be  independent  of  the  number  of  the  spectral 
bandi  filters  used  in  Fourier  plane.  In  other  words,  the  input  grain-noise 
is  considered  as  a  part  of  input  object,  which  can  not  be  eliminated  by 
temporally  partial  coherent  or  incoherent  illumination. 

4.3  Effect  of  Noise  at  the  Fourier  Plane 

We  shall  now  describe  the  effect  of  noise  at  the  Fourier  plane.  Sine 

photographic  plates  are  usually  used  to  synthesize  the  optical  spatial 

filters  for  signal  processing,  the  phase  fluctuations  and  granularity  of 

the  spatial  filters  have  significant  effect  on  noise  performance  of  the 

18 

optical  processor.  Although  the  setup-noise  is  primarily  due  to 
scattering  and  reflections  from  the  optical  components  of  the  optical 
processing  system,  however  we  assume  that  all  of  the  optical  components 
are  sufficiently  far  from  the  input  plane,  as  compared  to  the  focal 
depth  of  the  optical  system.  It  is  therefore  reasonable  to  assume  that 
all  the  setup-noise  sources  are  concentrated  in  the  Fourier  plane.  If 
the  optical  components  are  of  high  quality,  the  resulting  setup-noise 


sed  sii 


is  generally  very  small  as  compared  with  the  phase  and  the  grain-noise 

of  the  photographic  emulsion.  Thus,  the  phase  and  grain-noise  of  the 

nth  spectral  baud  filter  11  (0)  may  be  written  as, 

n 

Hn((0  =~[exp(i  K$((S))+  W)]  (*^y)  ?  <77> 

where  (8 )  and  G(0)  are  the  respective  random  phase  and  grain-noise 
distributions,  the  rectangular  function  represents  the  size  of  the  spatial 
filter  in  the  0-direction,  which  is  limited  by  the  spatial  frequency 
of  the  input  signal,  v.  Again,  for  simplicity,  we  would  use  a  one¬ 
dimensional  notation  for  the  noise  analysis. 

We  shall  now  analyze  the  effect  of  the  noise  performance,  due  to 
the  random  phase  <j>(6)  and  grain-noise,  respectively,  G(0)  at  the  output 
plane  of  the  proposed  white-light  optical  processor  in  the  following: 

4.3.1  Phase-Noise  _ 

Again  we  utilize  a  sinusoidal  grating,  as  input  signal,  superimposed 
with  a  phase  grating  at  the  input  plane.  The  resultant  amplitude  trans¬ 
mittance  is 

•fc(x, jj)  3—( (4  CoSZnuy)  exp (CzirVoX.)  ,  (78) 

The  corresponding  Fourier  spectrum  at  the  spatial  frequency  plane  can  be 
written  as 

S(oi-D0f\(5)-5-(a-i;.fx)p((3)+  >  (79) 

With  reference  to  Eq.  (39)  the  ensemble  average  of  the  output  intensity 
due  to  the  nth  spectral  band  filter  would  be 


Etwru-S 

'°° 

2 

+  X^((iz-i;5x)4  1- £( 

-  hea.  (  t-M*  (-jpr)  E  [e*Pl‘K  W') 

-4>((J2)3}  3  e*p{-CK  |j-C(<=i,-«ta)X+  C(VP0tf3] 

<  4.o*i4d2  4/Si  4(5*  4x.  . 


From  the  result  of  Eq.  (61)  Eq.  (80)  can  be  written  as 


E  fines'))  =AX*(H  ccs2nvs')j-  <**•*  {* 


-  4  K*(AAn)3  CoS2*W+  [ 


inhere 


-4  K,  (aX*)]}  6054*^'  > 


A**t  (  A^nt  A\*) 


-  . 

K,(4An)Aj  -^-expC-Mn/cOdA 


—  Vfn  exP  (“^n W*  )  ~  \jjri  e*P  (-  *jtn  fV/^) 


A*n 

Ki(M»)  — j  -jr  exp(-\i»/zd.)ix 


\m 


_  XE,g*P(-xXnW^)a^exp(-AM*>V*<l) 
^■■An  ^JUv 


-|rfE t(-WwM)*E«(-\u,$>>/*4,)]  , 


Similarly,  the  ensemble  average  of  I  (y * >  can  be  evaluated  in  the 
following: 

EU’o'))  -niH  ^(P0*rs<'5'-^x'>zs’((S'*vfx')] 

-  oo 

« i  s'  -  uf\ty  2  r(p ,+ vi  xj)] 

♦  K.Ki[e[<t>(W4(Pii]',E[4,(?0<*“(Pt)] 

+  E  [4>(|>»)4>((5,)]  -  E  [4> ((**)<♦> ((*03 

•»  kJ  e  [4>(Pi)«<>(po]  *■“■*  (  ^7-) 

*  «**  ( ifs -)  *“* 

*exp{-£  3^™  ((*r (**)$- <*  3^’(?rP*)J  } 

x  4(5, 4.^-1  tflAz  .  (8 

By  substituting  the  exponential  form  representation  of  Eq.  (53)  i.e., 

E[  d>  (Bj^)  d>  (62)  ]  =  °2  exP[_|S1~82|/dJ  in  Eq.  (82),  we  have, 


e  lie?, 


3  A\> 


+4«Vax,(i»«s*. vp  [x-(^-,xO 


+  Ky(AAn)c0SE7TV»^\  Kg  (&\n)C0$4*V  jj'  J 

+  K3 (^A«)(K4(AA.n^-Kj(AX«-}3 

*  (I-  cos v y' ) (\  +  cjos 2ii iJy')1, 
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where 


XjLn 
A  f  l 


rHn 


(83a) 


K4(aAn)4^  -i- exp ( f  Xi)/ 2.4.)4A= E i(\juW2^)~Ec(^^/2<^) 

Xa*  (83b) 

Ks(&x*)  4  2.K*  (AX,)  -  —^4 - Op(-Aj.Aj-D/2J.) 


-  exp(-x,ujy/2a)]  ; 


(83c) 


*UtA 


K6(M,)=K,(aA,)/2-  A 


i  V  ^Jtn  ^4ln 
-  exp(-X^fv/<il)]  . 


[«*P(-WvM) 


(83d) 


With  reference  to  the  definitions  of  Eqs.  (41)and  (42), the  output  SNR^, 
due  to  nth  spectral  band  filter,  under  a  broadband  temporally  partial 


coherent  illumination  is 


snr'”=  AJU  {( I*  cos^^ttV 

,  ccszvny  >  J^aX.)'  4  / 

{/in20'2(  |-f  coSZ*^)*  [K3(AA")(K4(a\,0-  K5  (<aO)/AXn 
+  (KSC^)-^Kl(AA.n)-f  -  Alb^L  )C0S2*V'4/ 

4[K6UAn)-K,(AXn) - - K3(AX»\) 


X-ju.  CXa* * ^Xn) 

X  [k4  (ax*)-  Kj  (AXn)j  /a  a  n  J  cos  4  *  y  g '  ]  ),/2 


(84) 


where  the  superscripts  (3)  represents  the  effect  due  to  phase-noise  at 
the  Fourier  plane.  Futhermore,  we  note  that 


(84a) 


K I  (AXu)  f  to\r 

tX*+Q 


__  I 


A.tn  A-M  ' 


^•v  K?  (AAn)/AA» 

AA„-»0 


A&»\  Am  2 


) 


(84b) 


JLi*n.  Kj(AA«)/dAn 
A\„-»0 


A  |Ui  2^. 


(84c) 


jU**  K$  (AAn)/aAn  — • 


AArH^O 


JL.4  ML 

Aju  2^ 


(84d) 


JU^n  |<j(AA*)/AAn 

4Xii'+0 


<81B) 


2 
A j* 


and 


JU-v  Ke(AX„)/«n  -  -j-(  ir-x)+  ■k6xp('X^iv/‘A)- (84f> 


2AjUx'  At*  «.  ' 

4X»i+0 

Thus,  the  SNRn  for  strict  temporally  coherent  illumination  would  be 


SNR«VW~t  =JU~  5N""V> 

CotxKMnX  AXm*0 


=  {(  l+OOS2Trv4')2-4nVl['Tr'Cc527T2)^  +  — s-  C0S4*»V)']}/ 

+  -5Tr«p(-w>»/**))«-s»^V((  J-  -  *z*  -J-r 

^“jU 

€*p  c- Am  f^M))ws4n^^']  y/z 


+  -±±L  ♦ 


I 


15T 


(85) 


Finally,  the  normalized  output  SNR,  due  to  nth  spectral  band  filter,  can 
be  shown  by  the  following  definition: 

. 

1 

mm  l 


5NR„  —  M  £  5NR*  (^)/SNRn  c^w)  ^  (86) 

(•<■1 

Since  conclusive  results  are  rather  difficult  to  deduce  from  Eqs.  (85) 
and  (86)  we  would  convert  these  equations  in  more  amenable  forms,  for 
which  are  suitable  for  numerical  computation.  Let  us  confine  the 


normalized  spatial  frequency  of  the  input  object  transparency  in  a 
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region,  that  the  following  condition  holds: 

.n  <  <s 

where  d  in  Che  correlation  distance,  f  Is  the  focal  length  of  transform 


lens,  X  is  the  wavelength  of  illumination,  and  v£  is  cutoff  spatial 
frequency  of  optical  system. 

We  note  that  the  correlation  length  d  of  the  phase  fluctuation 

can  be  assumed  equal  to  the  minimum  resolution  distance  of  the  optical 

processor.  Equation  (85)  is  essentially  representing  a  low  spatial 

frequency  signal,  since  the  ratio  of  d/  (fXv^)  is  generally  true  for  the 

proposed  white-light  optical  signal  processor.  Thus  we  shall  use  the 

condition  of  Eq.  (87)  [i.e. ,  fl  <  d/ (fXv  )]  to  describe  the  noise 

c 

performance  due  to  phase-noise  at  Fourier  plane.  By  the 

imposition  of  Eq.  (87)  Eq.  (86)  can  be  written  in  the  following  form, 

which  is  more  suitable  for  numerical  computation: 


x  (CoSzm>^m  +  OOStt«V^)]{\-$coS2in>y'm 

-  2  cos  4  ]'„)/  {  [ ( I  ♦  CoS  Z*V  fli,)2- 

x  (2cx»s2TTi;^  +  co54ny^l,)][(i  +  co527ryjl,) 
x  hm  Mh>M./^)  ,(2eoSW^ 


+  COS  4W„)  ^  ]1 

®  Xi*,AX«  -1  • 


From  the  description  of  Eq.  (88),  and  note  that  the  term  [£n(l  +  AX/X.  ]/AX 

£n  n 

is  a  monotonically  decreasing  function  of  AX^,  we  see  that,  the  broader 

-  (3) 

the  spectral  bandwidth  AX  ,  the  higher  the  SNR  can  be  obtained. 

n  n 

Since  AX^  is  inversely  proportional  to  the  spatial  frequency  vq  of  the 

-  (3) 

phase  grating,  a  higher  SNR^  could  also  be  obtained  by  using  a  lower 
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v  .  However  the  number  of  spectral  band  filters  also  reduces,  for 

a  lower  vq.  We  should  emphasize  that,  the  larger  the  number  of 

- (-1} 

upeetrul  hand  flltoru  unod ,  the  higher  the  overall  out|)ut  SNK  can 


be  obtained,  although  SNR 


(3) 


decreases  somewhat  Cor  higher  v  .  T.n 

o 


otherwords,  the  filtered  signals,  from  each  of  the  spectral  band  filters, 
would  be  incoherently  superimposed  at  the  output  plane.  Thus  the  effect 
due  to  the  phase-noi6e  from  the  spectral  filters  Hn(0)  would  be  sub¬ 
stantially  supressed  due  to  the  incoherent  addition.  Thus,  the  overall 
“—(3) 

output  SNR  ,  due  to  phase-noise  at  the  Fourier  plane,  can  be  written 

as  ,  N 


n- 1 

We  note  that,  by  replacing  AA^  by  AA  in  Eq.  (89),  we  have 


(89) 


-  v  - 

5MRWr*  SNR* 


(90) 


where  AA  ■  NAAn>  the  entire  spectral  band  of  the  white-light  source. 

We  also  note  that  this  results  is  quite  consistent  with  the  experimental 

results  for  broadband  deblurring  with  a  white-light  source^. 

In  order  to  get  some  more  feelings  on  the  effects  due  to  signal  spatial 
(3) 


frequency  on  the  SNR 


,  we  let  AAfl  go  to  infinity,  i.e., 


'  AXn  *  '  * 


°°  OA* 

Thus  Eq.  (88)  reduces  to  the  following  form 


=  |  (91) 


_  .  M 

SNR*5  ==  {(t4Co$2ituy'MXh?cc$2iJljy/m+2C0$4niJy,„))/ 


(92) 


^..  .  V-A  i  ■  %  V- 
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-  (3) 

Thus,  we  see  chat,  for  relatively  low  ft,  the  SNR^  increases  as  the 

spatial  frequency  of  the  signal  increases.  In  other  words,  a  significant 
-  (3) 

improvement  of  SNR^  can  be  obtained  for  higher  spatial  frequency  signal. 
However,  the  higher  the  spatial  frequency  requires  a  larger  size  of  spectral 
band  iilter  [i.e.,  AH^fB)].  For  temporally  partial  coherent  illumination, 
the  noise  at  the  output  plane  would  also  be  increased,  compared  with  a 
lower  spectral  signal.  However,  the  situation  is  quite  different  for 
multi-spectral  band  filterings,  the  overall  output  noise  would  be  sub¬ 
stantially  supressed,  due  to  the  incoherent  addition  of  different  spectral 

band  filtered  signals. 

-  (3) 

Several  plots  of  SNR^  as  a  function  of  spectral  bandwidth  AAn> 

for  various  values  of  normalized  input  spatial  frequency  ft,  are  given 

in  Fig.  2.0,  The  plot  for  low  frequency  case  (i.e.,  ft  =  0.1)  was  evaluated 

by  using  the  Eq.  (92)  while  for  ft  =  0.4  the  plot  for  ft  *  0.4  is  obtained 

-  (3) 

by  using  Eq.  (88).  From  this  figure  we  conclude  that  the  SNR^  ,  due  to 

phase-noise  at  the  Fourier  plane  can  be  noticeably  improved  by  using  a 

broadband  temporally  partial  coherent  illumination,  which  is  similar  to 

the  effect  due  to  the  phase-noise  at  the  input  plane,  as  shown  in  Fig.  4. 

-  (3) 

The  higher  the  spatial  frequency  of  the  input  signal,  SNR^  can  be 

improved  to  several  tens  of  orders,  as  compared  with  low  spatial  frequency 

-  (3) 

case.  We  note  that  the  improvement  of  SNR^  is  several  orders  higher 

than  SNR  duu  to  the  phase-noise  at  input  plane,  for  higher  spatial 

-  (3) 

frequency  case.  Similar  to  the  input  phase-noise  effect,  the  SNR 

n 

can  also  be  improved  for  shorter  lower  wavelength  limit  .  The 

-  (3) 

variation  of  SNR  as  a  function  of  spectral  bandwidth  AA  ,  for 
n  n 

various  values  of  A  ,  are  plotted  in  Fig.  21.  From  these  results,  we 

-  (3) 

see  that  the  output  SNR  is  approximately  exponential  increasing  as  a 

n 

function  of  AA^.  As  similar  to  the  effect  due  to  phase-noise  input 
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plane,  the  improvement  is  somewhat  higher  for  shorter  wavelength  temporally 
partial  coherent  illumination. 

4.3.2  Grain-Noise 

In  the  evaluation  of  the  noise  performance  due  to  grain-noise  at  the 
Fourier  plane,  we  shall  handle  the  problem  somewhat  differently.  Since 
the  scale  of  the  signal  spectrum  is  proportional  to  ‘the  illuminating 
wavelength  of  the  light  source,  the  diffraction  due  to  each  of  the 
spectral  band  filters  behave  quite  differently  with  different  wave¬ 
length  illumination.  The  grain-noise  at  the  filter  plane  can  be 
substantially  smoothed  out  with  a  broader  band  temporally  partial 
coherent  illumination.  Thus,  we  would  expect  a  higher  output  signal- 
to-noise  ratio  for  a  broadband  illumination  (e.g.,  white-light).  We 
shall  now  provide  a  quantitative  analysis  of  the  effect  due  to  grain- 
noise  at  the  Fourier  plane.  Here  we  utilize  a  narrow  spectral  band 
filter  with  an  additive  film-grain  noise,  as  described  in  the  following: 

.  (93> 


By  substituting  this  filter  function  into  Eq.  (39),  and  employing  the  same 
analysis  as  we  did  for  the  phase-noise  case,  the  mean  value  of  the 

output  intensity  distribution  can  be  determined  by  the  following  integration: 

oo 

E  (In (g')j  =  E  (  JJ  JJ[J  S A.,  (3i)S‘<c<*-M>f\  p«) 

*  (i+CiCP'))  [i+dcpo] 

%  d(3| 


(94) 


With  reference  to  Eqs.  (81),  (49),  and  (50),  the  above  equation 


can  be  written  as: 

oo 

Ajtn 

4  4-r((5j)S(fJ,tufA.)4  -bff(i,-^iA)£-((st-ufx) 

+  -l-r(Pi*VfX_)5(f!j4yfX)4-i-S'((?,4jjfx) 

*  [ |  +  2CAtCA1t£J.p(2.3D>  I ] exp [-1(1, -p, I /«.] 

tCA}e*P[-i^-((3rp03,]a(J.a(JzilX  >  (,5) 

which  is  equal  to 


£  ®*AAn(|4CA)2(/4C052»r^j/)2+  [exp(llD) 

-l][3AX«  +  K0(a\h)  CoS 2nvy' 

W*i.K.  +Ko(AAn)C0S4Tryj')]  f 

K0(AXn)  ^•exp(-Ajln5-^/2a')[*‘exP(~^r,i'i;/2A')J  ,  (g6a 


(96) 


and 


Ko(AAn)  *exp(-AJtw5v/a.j)(|-exp(-^Anfi>/^)]  . 


(96b 


For  the  ensemble  average  of  the  square  irradiance,  we  have 

E  (if(g')]  -EtH*  flfl  t5CP*)+iff(Pr»fA,)4^S(pl+vJx,)] 

A**  -«• 

*  [«'<(J*)  +  -z5'C(3i- vfA,^4-4-5"(fJe+  x>fX,)) 
«(5-(fi,)*T5'((4-"^)4x*'(P‘  +  >'JX0] 


x  i  s  (fa  -  vf  a2)  4  2“  x0] 

+  $,((*,)  G> ( (f,)  4  G,((5.)  <*(*) 4 C1  ( (*0 °I  <  (*•  )♦  fyPOWO 

4  Cj(/S.)C-)C?i)]  expi-c  x^(P«-(U)^-<-^j((V 


By  substituting  the  mean  value  and  the  autocorrelation  function  of  the 
film  granularity  (same  as  Eqs.  (49)  and  (50)  into  Eq.  (97)  the  ensemble 
of  the  square  irradiance  would  be 

E  (ii  c p]  =  AA*  ( I*  4 CA 1 6 c\  )( u co&2n»y‘ ')* 

+  2Cl(exp(2-3D)-  l)(|  +  c«S2irys')J 

*^,tf£&;K*(“,0eo5*,ni 

2 

4  K7(A*n)  (98) 

jzvx*&  0  J  > 

where 

K7  (  AXn)  JL  exp  (A An  f  ^/2  a)  4  exp  (-AAn^/2^) 

-$4  exp(X^ 

-  z  exp  [-  (^Xjtn  4  &Xrx)  /  Za»]  m  (98a) 

From  the  results  of  Eqs.  (96)  and  (98),  the  output  SNR,  due  to  nth 
spectral  band  filter,  for  a  temporally  partial  coherent  illumination 
can  be  obtained  such  as 

5N^>Cp={(l*CA),(|*toS2WVp%-i-ci(e><PU.30)-iJ 

x  f  3  ♦  —  ■■■  (8Ko(&*n)coS2ir»'fi  Ko(AXn) 

^  JV AAn x  V 
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xcos4ttvY  )]}  /  {Ca  (i+  cosaitv^)2  [e*p(2.3i>)- 1  ] 

*  C«  4  K.  M  ccszn^'f 

K7  («*•*)- *.'  (^.))eoS«iriij']}'/'1  <”> 


where  the  superscript  (4)  denotes  the  effect  due  to  grain-noise  at  the 
Fourier  plane.  It  is  also  interesting  to  note  that,  by  letting  AA^ 
equals  to  zero,  we  obtain  the  following  output  SNR  for  strictly  temporal 
coherent  illumination,  i.e., 

SNft')(PcoWJt=t(|4CA)l(l*co32,ri;^:)il+^ 

X  [exp(2.3D)-  I  ]  [3+4exp(-Ai«,f>'/Z«.)C0S2»ri;g' 

+  exp  (-  “54"  »  jj'  J )  /{<*  ( I + Cttsm»pl 

x  (exp (?-3D>  I  ]  [4  *  6 exp  (-  A**  f  u/20.) CCS 2nW 
- -1-cxp  frr  AXn»0  (100) 


From  Eqs.  (99)  and  (100),  the  normalized  SNR,  due  to  nth  spectral  band 
filter,  can  be  obtained,  such  as 

_ .  M 

SNRn  “ITZ  {0’ca)  ll'*<:*S2nVjl,)%-l<j|(e,tP(2’3P.) 

Mm| 

-  1 31 3  ♦  JZZZt  t 8  Ko  COS  ZlT  V  ^  +  kJ  (Mn) 

* co$4irv«^)j  [4  +  &exP(-xJz*5-^/2a)c0S2nv^ 


f  ^  4  filTn  *  4 

jrW  AAn  w 

*  [(|  +  CA)2(l  +  C052Tri)^ln)2+  -^“CA  [€*P(2*3D) 

-I]  [3  +  2  exp(-XJUxf'^/20-)coS2TTX)^ 

+  exp(-XJln^/a-)cos4n^^]}  ?  (101) 

where  we  note  that 

Xi«.  K.(4>-n)/4An=>|^exp(.Airtfi'/2<».)  ,  (1014) 

AAn->0 

JUUn.  Ko(*x*)M*  €Xp(-X^f»/a.)  ,  (101b) 

AA»-*0 

JLtm  K7  (AAn)/ &An  *  0  .  (101c) 

aA*-»  o 

We  stress  that,  for  all  practical  cases,  the  mean  value  of  the 

2 

grain-noise  is  much  smaller  than  unity,  (i.e,  C  «  C  ).  Thus  Eq.  (101) 

A  A 

can  be  reduced  to  a  simpler  form,  such  as 

5NR^=-L-£  {4t6exp(-A >jUtW24.)c0S2irUjj^ 

-TcxP(-v*v/^ws«"vtf~l/l«+^n 

x  Ko(AAn)  COS2«V^  +  (  j£rn  *1  ^ 

-  Ko(4Xn)co54Tt^^  )J  (102) 

For  low  frequency  signal  [i.e.,  Eq .  (87)  ,  Eq .  (102)  can  be  further  reduced 

SNR1*’  =  <103) 
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-  (4) 

From  the  results  of  Eqs.(102)- and  (103),  we  see  that  the  SNR^  Increases 

as  the  spectral  bandwidth  AA^  increases.  This  means  that  a  higher 
- (4) 

Sh’R^  is  achievable  with  a  broader  band  temporally  partial  coherent 

-  (4) 

illumination,  as  predicted.  Several  curves  of  the  SNR^  as  a  function 

of  spectral  bandwidth  AA^,  for  various  normalized  signal  frequencies 

(0  a  v/vc),are  plotted  in  Fig.  22.  From  this  figure,  again  we  see  that, 
-  (4) 

the  SNR^  increases  exponentially  as  a  function  of  spectral  baridwidth 

AAfl  of  the  filter,  for  higher  spatial  frequency  of  the  input  object. 

On  the  other  hand,  if  the  spatial  frequency  of  the  input  signal  is  low 
-  (4) 

(e.g.,  v  -  0),  the  SNR^  would  be  independent  of  the  spectral  band¬ 
width  AA  .  That  is,  the  output  SNR  would  not  be  improved  by  using  a 
n 

broader  band  temporally  partial  coherentillumination.  We  stress  that  this 

behavior  is  quite  different  from  the  effect  due  to  phase-noise  at  the 

-  (3) 

Fourier  plane  for  which  a  higher  SNRn  can  always  be  obtained  by 
increasing  AAn,  even  for  zero  spatial  frequency  at  the  input  object 
(i.e.,  0  =  0),  that  is, 

(104) 

It  is  also  worthwhile  to  show  the  asymptotic  results  of  SNR"r/  ,  due  to 
phase-noise  at  Fourier  plane  for  low  frequency  signal,  i.e.. 


From  these  asymptotic  results  of  Eqs.  (105) to  (109), we  confirm  our 

conclusion  stated  previously;  nam  ly  the  SNR^  '  can  always  be  improved 

for  phase-noise  at  the  Fourier  plane  for  a  broader  band  temporally 

partial  coherent  illumination,  while  for  grain-noise  at  the  Fourier 

-  (3) 

plane,  the  SNR^  can  not  be  improved  for  v  *  0.  However,  for  high 

frequency  case  (e.g.,  Q  *  0.4),  the  SNR^  increases  quite  noticeable 

as  AX  increases  as  shown  in  Fig.  22. 
n  b 

-  (3) 

Although  SNR^  ,  due  to  grain-noise  at  input  plane,  is  monotonically 

decreasing  function  of  the  mean  value  CA  of  grain-noise,  as  shown  in 

-  (4) 

Fig.  8,  however  the  SNRn  is  independent  of  C^,  i.e.. 


- - (A) 

5NRn 


1  +  Ca 


for 


(108) 


which  is  obtained  from  Eq.  (67). 

We  shall  again  note  that,  the  overall  normalized  output  signal-to- 
- (4) 

noise  ratio  SNR  can  be  obtained  by 

SNR1*0  —  —  y  5noc<) 


irL  SNR"  , 


where  N  is  the  total  number  of  the  spectral  band  filters,  at  the  Fourier 

- (4) 

plane.  Thus  a  higher  SNR  can  always  be  obtained  for  larger  N,  particularly 

for  higher  input  spatial  frequency  case  (i.e.,  0  >>  0).  We  further  note 
-  (4) 

that,  although  SNR  is  independent  of  AX  as  shown  in  Fig.  22.  however 
n  n 

- (4) 

SNR  would  have  significant  improvement  since  the  output  noise  distribu¬ 
tions  (i.e.,  grain-noise  at  Fourier  plane)  due  to  each  spectral  wavelength 
X,  are  mutually  incoherent.  In  other  words,  the  signals  diffracted  from 
the  spectral  band  filter,  would  be  precisely  superimposed  at  the 
output  image  plane  to  form  a  higher  intensity  signal,  while  the  noise 
diffracted  from  the  filter  plane  would  be  randomly  added  to  form  a 
smoother  noise  background.  Thus,  a  higher  overall  output  .  .«R  can  always 


be  obtained  due  to  grain-noise,  as  well  as  phase  noise,  at  the 
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Fourier  plane.  Nevertheless,  the  number  of  the  spectral  band  filters  is 
generally  limited  by  the  spatial  frequency  of  the  phase  grating. 
However,  the  higher  vq,  the  larger  the  size  of  the  achromatic  transform 
lenses  are  required.  But  larger  achromatic  transform  lens  Is  generally 


expensive  and  more  elaborate  to  construct  in  practice. 


4 . 4  Summary 


In  summary,  we  would  point  out  that,  the  effects  due  to  temporal 
coherence  on  the  noise  performance  of  a  white-light  optical  signal 
processor  have  been  quantitatively  evaluated  from  partial  coherence 
theory  stand  point.  Both  the  phase-noise  and  the  grain-noise  at  the 
input  and  Fourier  planes  are  treated  in  the  evaluation.  We  have  shown 
that  these  noises  have  significant  effects  on  the  system  noise  performances 
with  temporally  partial  coherent  illumination,  as  summarized  in  the 
following: 

1)  There  are  significant  effects  on  output  noise  reduction,  due  to 
phase-noise  at  the  input  and  Fourier  planes  and  grain-noise  at 
the  Fourier  plane,  by  a  broadband  temporally  partial  coherent 
illumination. 

2)  There  is  however  no  significant  effect  on  noise  reduction  due 
to  grain-noise  at  the  input  plane  by  a  broad-band  temporally 
partial  coherent  illumination. 

3)  Except  due  to  the  grain-noise  at  the  input  plane,  the  output 

normalized  signal-co-noise  ratio  SNR^  are  generally  monotonicallv 
increasing  functions  of  the  spectral  bandwidth  of  nth 

spatial  filter. 

4)  The  output  normalized  signal-to-noise  ratio  SNRn  increases  (except 
due  to  grain-noise  at  input  plane)  as  the  input  spatial  frequency 
increases.  However,  increases  the  input  spatial  frequency  would 


also  increase  the  spectral  bandwidth  of  the  filter.  Thus 


it  reduces  the  number  of  the  processing  spectral  band 
filters  at  the  Fourier  plane.  Tor  a  fixed 

5)  The  normalized  slgnal-to-noisc  ratio  SNR  also  Increases 

n 

(except  due  to  grain-noise  at  the  input  plane)  as  the 
lower  wavelength  limit  (i.e.,  A  )  of  the  spectral  band 
filter  decreases.  In  otherwords,  the  SNRn  improves  for 

shorter  wavelength  temporally  partial  coherent  illumination. 

-  (2) 

6)  The  normalized  signal-to-noise  ratio  SNR^  is  a  mono- 
tonically  decreasing  function  of  the  mean  value  CA  of  grain- 

-  (4) 

noise  at  the  input  plane.  But  the  SNR  ,  due  to  the 

n 

grain-noise  at  Fourier  plane,  is  independent  of  C^. 

7)  The  overall  normalized  signal-to-noise  ratio  SNR  is  equal  to 
the  sum  of  SNRn»  except  for  the  grain-noise  at  the  Input 
plane  which  is  independent  of  n.  SNR  Improves  as  N  number 
of  spectral  band  filters  increases,  as  shown  in  Fig.  16. 

Finally,  we  conclude  that,  the  proposed  optical  signal  processor 
i3  capable  of  improving  the  noise  performance,  with  a  broadband  temporally 
incoherent  (i.e.,  white-light)  source.  Except  the  grain-noise  at  the 
input  plane,  the  technique  is  capable  of  suppressing  the  noise  effects 
at  the  input  and  Fourier  planes,  as  an  incoherent  processor,  although  the 
processing  is  carried  out  by  a  temporally  partial  coherence  mode. 

V.  Polychromatic  Matched  Filtering 
5.1  Polychromatic  Correlation  Detection 

We  shall  demonstrate  a  grating  based  polychromatic  coherent  signal 
correlation  detection.  We  shall  show  that  this  polychromatic  correlation 
detection  technique  is  simple  and  flexible  to  operate  and  this  technique 
offers  the  advantage  of  multicolor  signal  detections,  which  is  a  step 
closer  to  the  actual  human  visual  recognition. 
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Let  us  insert  a  color  signal  transparency  in  contact  with  a  sinu¬ 


soidal  diffraction  grating  at  the  input  plane  of  a  polychromatic 

coherent  optical  professor ,  as  shown  in  Fig.  23.  The  complex 

light  amplitude  distribution  behind  the  input  plane  can  be  described  as 


s(x,y;X)T(x,y)  =  [sr(x,y;Xr)  +  sg(x,y;Xg)  +  sb(x,y ;Xb) ]T(x,y) ,  (109) 

where  s  ,  s  ,  and  s,  are  the  red,  green,  and  blue  amplitude  transmittances 
r  g  d 

of  the  color  signal  transparency,  X^,  X  ,  and  Xb  are  the  red,  green,  and 
blue  wavelengths  of  the  coherent  sources,  and  T(x,y)  is  a  sinusoidal 


grating  which  can  be  described  as, 
T(x,y)  -  K[l  +  cos(pox)] , 


where  K  Is  a  proportionality  constant,  and  is  (.he  angular  spatial 
frequency  of  the  grating. 

Since  the  input  plane  is  illuminated  by  a  collimated  polychromatic 
coherent  plane  waves  of  red,  green,  and  blue  wavelengths,  the  complex 
light  distributions  at  the  back  focal  plane  of  the  achromatic  transform 
lens  are 

fX 

S(ct,6  ;X)  -  Sr(«,8;Xr)  +  Sg(a,8;Xg)  +  Sb(a,8;Xb)  +  Sr(a  +  ~  pQ,e) 


+  Sg(a  +  Po,8)  +  Sb(a  ±  ~  , 


(110) 


where  we  have  ignored  the  proportionality  constant  for  simplicity, 

(a, 8)  denotes  the  spatial  coordinate  system  of  the  Fourier  plane,  f  is 

the  focal  length  of  the  achromatic  transform  lens,  and  S  (a, 8),  S  (a, 6), 

and  S, (a, 8)  are  the  corresponding  color  Fourier  spectra  of  s  (x,y), 

D  r 

s  (x,y)  and  s,  (x,y).  From  this  equation,  we  see  that  the  red,  green, 

8  D  fX  fX 

and  blue  color  spectra  are  centered  at  a-+  -r—  p  ,  a=+  p  ,  and 

ZTT  O  Zn  O 

^  ^  b 

a*+  PQ  locations  along  the  o  axis,  respectively.  Since,  in  practice, 
the  all  color  signal  transparencies  are  spatial  frequency  limited,  it 
is  possible  to  calculate  the  appropriate  spatial  carrier's  frequency 
Pq  of  the  grating  so  that  the  primary  color  Fourier  spectra  would  be 
spatially  separated  in  the  Fourier  plane. 
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In  color  signal  detection,  we  assume  that  a  set  of  color  sensitive 


The  corresponding  amplitude  ‘ 

i 

fA 

axo  +  2 IT  Po*6)1* 

fA  ^ 

aXo  +*g(a-  2^P0*S)]* 

f>,  ' 

axo  +  V°“  2T  PQ*e)J  (111)  j 

frequency,  Ks  are  the  appropriate 


(112) 

are  the  primary  color  detecting  signal  spectra. 

If  we  properly  insert  these  color  sensitive  matched  spatial 
filters  in  the  Fourier  plane  p^  of  the  polychromatic  coherent  processor, 
as  illustrated  in  Fig. 24*  the  complex  light  distribution  at  the  output 
plane  p^  of  the  processor  would  be 

g(x,y)  *  E[sn(x,y)exp(ipox) 

+  sn(x,y)  exp(ipQx)  *  sn(x  -  xQ,y)  exp(ipQx) 

+  sn(x,y)  exp(ipQx)  *  sn(“x  +  xo>y)exp(ipox) ] ,  for  n»r,g,  and  b,  (113) 

where  *  denotes  the  convolution  operation,  and  we  have  ignored  the 
proportionality  constants  for  convenience. 

From  the  above  equation,  we  see  that'  the  first  summation  terms 


matched  spatial  filters  are  constructed, 
transmlttances  are 


fA 

fX 

Hr(o,8) 

-Kl  + 

K| 

S  (a-  ~ 
r  2  it 

Pa.8) 

Ic°s[^r 

fX 

fX 

Hg(o,6) 

-Kl  + 

K| 

V"  2^ 

P0,8) 

|cos[^ 

fA,_ 

fx. 

^(o.B) 

-Kl  + 

k| 

lSb‘“-  IT 

Po’8) 

I«*l5k 

where  x 

o 

Is  an 

arbitrary  carrier 

spatial 

proportionality  constants. 

and 

sr 

(a, 6)  ' 

- 

|Sr(a,B)  | exp[ i  q> 

r(a,3)]. 

ss 

(a,0)  ’ 

- 

|Sg(o,8)  |exp[i  <j> 

g (a,6)3. 

and  Sb(a,B)  -  js^a.B)  |exp[i  $b(a,6)J, 


are  the  zero-order  terms  which  are  diffracted  in  the  neighborhood 
of  the  optical  axis  at  the  output  plane,  and  the  second  and  third 


Multicolor  complex  spatial  filtering.  H 
are  red,  green,  and  blue  color  sensitive 
spatial  filters. 
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stunm.it  Jon  terms  are  the  primary  color  convolution  and  correlation  terms, 

which  are  diffracted  in  the  neighborhood  of  (-x  ,0)  and  (x  ,0)  at  the 

o  o 

output  plane  of  the  polychromatic  coherent  processor,  as  sketched  in 
Fig.  25.  From  this  result,  we  see  that  this  proposed  polychromatic  pro¬ 
cessing  technique  is  capable  of  performing  color  signal  correlation 
detection  and  the  technique  is  rather  simple  and  flexible  to  operate. 

One  should  also  note  that  this  polychromatic  process  technique  is 
also  capable  of  performing  multicolor  signals  detections,  if  the 
conjugate  orders  of  color  Fourier  spectra  are  utilized,  as  depicted 
in  Fig.  26.  Where  the  subscripts  of  1  and  2  represent  two  different 
sets  of  matched  spatial  filters  for  color  signals  no.  1  and  no.  2, 
respectively.  We  would  also  note  that,  if  each  of  the  color  signals  to  be 
detected  contain  only  one  primary  color,  then  it  is  possible  to 
synthesize  a  total  of  six  matched  color  sensitive  spatial  filters 
for  simultaneous  signal  detection.  Needless  to  say  that,  by 

19 

utilizing  a  multigrating  technique,  as  described  in  a  previous  paper  , 
this  technique  is  capable  of  performing  N  number  of  color  signal  corre¬ 
lation  detections,  with  relative1 ease  of  physical  operation,  as  illus¬ 
trated  in  Fig.  27. 

5.2  Color  Sensitive  Spatial  Filter 

We  shall  now  describe  a  technique  of  generating  a  set  of  spatially 
multiplexed  color  sensitive  matched  filters  for  color  signal  correlation 
detection,  as  illustrated  in  Fig.  28.  This  set  of  color  sensitive  matched 
filters  can  be  generated  sequentially  with  respect  to  the  color  wavelengths 
of  the  coherent  sources.  For  example,  the  red  color  sensitive  matched 
filter  is  generated  when  the  shutter  of  red  coherent  source  is  open  and  the 
recording  aperture  of  the  encoding  mask  MS  is  located  at  the  position 
where  the  red  color  Fourier  spectrum  is  exposed.  And,  similarly  for  the 


Figure  28.  Synthesis  of  spatially  multiplexed  color  sensitive 
matched  filters.  BS;  beam  splitters,  M;  mirror,  SF; 
spatial  filter,  MS;  recording  mask,  H;  holographic 
plane,  AP;  recording  aperture,  L;  achromatic  transform 


generation  of  the  green  color  sensitive  matched  filter  and  the  blue 
color  sensitive  matched  filter,  with  the  same  color  detecting  signal 


transparency.  In  this  manner,  we  see  that,  a  set  of  color  signal 

matched  filters  can  be  sequentially  synthesized,  without  physically 

overlapping  in  the  Fourier  plane.  Similarly,  if  one  takes  advantage 

of  the  c'njugate  orders  of  the  Fourier  spectra,  another  set  of  the 

color  sensitive  matched  filters  for  a  different  color  detecting 

signal  can  be  generated,  in  the  lower  half  of  the  spatial  frequency 

20  21 

plane.  Unlike  the  wavelenght-multiplexing  techniques  ’  ,  higher 

diffraction  efficient  color  sensitive  matched  filters  can  be  synthe¬ 
sized.  There  in  another  advantage  of  this  technique,  a  large  number 
of  multiple  color  sensitive  matched  filters  can  be  generated  as 
illustrated  in  Fig.  27.  This  large  capacity  multi-signal  matched 
filtering  program  is  currently  under  investigation. 

VI .  Experimental  Demonstrations 
6.1  Character  Recognition 

We  snail  illustrate  an  experimental  result  obtained  with  the 
incoherent  signal  detection.  Figure  29(a)  shows  a  set  of  alphanumer ics 
as  an  input  object  transparency.  Figure  29(b)  represents  the  character 
G  to  be  detected,  for  which  a  complex  spatial  filter  was  constructed 
with  the  interferometric  technique.  Figure  29(c)  shows  the  correlation 
detection  obtained  with  the  incoherent  correlation  detection  technique, 
where  a  narrow  band  mercury  arc  lamp  was  used  as  the  incoherent  light 
source.  Figure  29(d)  shows  the  electronic  scanned  of  the  correlation 
peak  obtained  from  the  result  of  Fig.  29(c).  For  comparison  we  also 
provide  the  correlation  detection  peak  obtained  with  coherent  tech¬ 
nique  as  shown  in  Fig.  30(a).  Fig.  30(b)  shows  the  corresponding 
electronic  scanned  result  of  Fig.  30(a).  With  reference  to  the 


>  X 


O  (N 


Character  recognition  with  incoherent  lif»ht 

(a)  Input  alphanumeric; 

(b)  Alphabet  G  to  be  detected; 

(c)  Correlation  detection  of  G; 

(d)  Electronic  scanned  correlation  peak. 


results  obtained  with  the  incoherent  source  as  shown  in  Fig.  29(d), 
wc  note  that  higher  signal-to-noise  correlation  peak  can  be  achieved 
with  the  incoherent  correlation  technique. 

6. 2  Detection  of  Pseudorandom  (PN)  Code 

The  basic  idea  of  using  pseudorandom  (PN)  code  for  the  spread 
spectrum  coranunication  is  to  encode  the  informational  signal  (either 
with  direct  sequence  modulation  or  frequency  hopping  or  the  combination 
of  both),  so  that  the  signal  is  capable  of  anti-jamming  or  interference 
during  the  transmission.  Since  the  spread  spectrum  communication 
is  a  large  time-bandwidth  signal,  the  optical  correlator  is  particularly 
suitable  to  decode  the  broad  spectrum  signal.  We  shall  have  to  demon¬ 
strate  an  experimental  result  that  an  incoherent  optical  correlator  is 
suitable  for  such  applications.  Figure  31(a)  shows  a  pseudorandom 
(PN)  code  as  an  input  object  transparency.  Fig.  31(b)  represents  a 
code  work  that  we  would  wish  to  detect.  A  complex  spatial  filter  of 
this  code  word  was  constructed  with  the  interferometric  technique.. 
Figure  31(c)  shows  the  corresponding  correlation  peak  obtained  with 
the  incoherent  processing  technique.  The  light  source  utilized  was 
a  75W  mercury  arc  lamp  with  a  green  interference  filter  wavelength 
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5461A  and  100A  bandwidth.  Figure  31(d)  shows  the  profile  of  the 
corresponding  correlation  peak  as  obtained  by  electronic  densitometry 
scanned. 

In  comparison,  we  also  provide  the  correlation  peak  obtained  with 
the  coherent  processing  technique,  as  shown  in  Fig.  32(a)  Figure  32(b) 
shows  the  electronic  scanned  correlation  profile  obtained  with  the 
coherent  method.  By  comparing  the  results  obtained  with  the  incoherent 
processing  technique  of  Fig.  31,  we  note  that  the  result  obtained  by 
the  incoherent  technique  offers  a  high  correlated  signal-to-noise  ratio. 


Detection  of  pseudorandom 
(a)  A  pseudorandom  code. 
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P  uedorandom  code  detection  with  coherence  source. 

(a)  Correlation  detection  with  coherence  source. 

(b)  An  electronic  scanned  correlation  peak  of  (a) 
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We  shall  now  illustrate  that  the  proposed  optical  correlator  is 
suitable  for  color  signal  correlation  detection.  The  significance  of 
color  signal  detection  is  a  step  closer  to  the  actual  human  recognition 
process.  We  note  that  the  proposed  grating  based  correlation  technique 
could  either  utilize  coherent  or  incoherent  source.  For  simplicity, 
our  experimental  demonstration  would  utilize  a  polychromatic  coherent 
source.  In  principle,  the  correlation  detection  can  be  utilized  by 
the  entire  spectral  band  of  a  white-light  source,  where  a  research 
program  in  under  investigation. 

Strictly  speaking,  coherent  light  sources  emitting  red,  green, 
and  blue  color  wavelengths  should  be  used  for  the  color  signal 
correlation  detection.  For  simplicity,  we  shall  use  only  two  of  the 
primary  colors  of  light  sources  in  our  experiments.  In  the  first 
experimental  demonstration,  a  set  of  color  English  alphabets,  as  shown 
in  Fig.  33(a),  is  used  as  input  color  object  transparency  whereas  the 
upper  row  of  these  alphabets  is  in  green,  middle  row  is  in  yellow,  and 
lower  row  is  in  red.  A  red  and  a  green  color  sensitive  matched  spatial 
filters  for  detecting  the  color  alphabet  A  are  generated  in  the  Fourier 
plane,  as  described  in  previous  sections.  A  HeNe  laser  emitting  a 
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6238A  red  light  on  an  argon  laser  emitting  a  5145A  green  light  are 
used  in  the  proposed  polychromatic  optical  processor,  as  shown  in 
Fig.  23.  Figure  33(b)  shows  the  color  signal  correlation  detection 
peaks  of  the  color  alphabet  A.  In  this  figure,  the  left  hand  corner 
correlation  peak  is  in  green,  the  middle  correlation  peak  is  in  yellow, 
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and  the  right  hand  corner  correlation  peak  is  in  red.  From  this  result 
we  see  that,  the  color  sensitive  matched  spatial  filters  possess  a 
high  degree  of  color  correlation  detection  selectivity. 

In  the  second  experimental  demonstration,  a  continuous  tone 
color  aerial  photographic  image  transparency,  as  shown  in  Fig.  34(a) 
is  used  as  input  object.  The  color  of  the  ground  terrain  is  generally 
yellowish  with  scattered  spots  of  light  green.  The  forest  and  some 
scattered  trees  apperas  to  be  redish  brown  and  yellow.  The  lakes  and 
ocean  are  in  dark  green.  Figure  34(b)  shows  a  section  of  terrain 
is  used  as  a  color  detecting  signal.  A  red  and  a  green  color  sensitive 
matched  filters  for  detecting  the  color  terrain  of  Fig.  34(b),  are 
generated.  Again  a  HeNe  laser  for  red  light  and  an  argon  laser  for 
green  light  are  used  for  the  color  signal  correlation  detection. 

Figure  34(c)  shows  the  correlation  detection  peak.  In  reality  this 
correlation  peak  is  in  yellow  color. 

For  our  third  experimental  demonstration,  a  continuous  tone  color 
image  transparency  of  a  campus  street  scene  in  The  Pennsylvania  State 
University,  as  shown  i-  Fig.  35(a),  is  used  in  the  input  object. 

Figure  35(b)  shows  the  color  objects  (i.e. ,  a  blue  van  and  a  red  stop 
sign)  used  as  color  detecting  signals.  Since  the  stop  sign  is 
primarily  in  red  and  the  van  is  in  blue,  therefore,  a  red  color  sensitive 
matched  filter  for  the  stop  sign  and  a  blue  color  matched  filter  for 
the  van  are  independently  constructed  in  spatial  frequency  plane.  In 
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this  experiment,  a  HeNe  laser  emitting  a  6328A  red  light  and  an  argon 
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laser  emitting  a  4880A  blue  light  are  used  for  the  color  signal 
correlation  detections.  Figure  35(c)  shows  the  color  signal  correlation 
detection  result.  The  left  hand  correlation  peak  is  in  blue,  which 


(c)  The  corresponding  correlation  peak  is  primaril 
yellow. 


The  corresponding  correlation  detections 
the  correlation  peak  for  the  van  is  blue 
and  the  correlation  peak  for  stop  sign  i 
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represents  the  correlation  detection  of  the  blue  van,  whereas  the  right 
hand  correlation  peak  is  in  red,  which  represents  the  correlation 
detection  of  the  red  stop  sign. 

In  concluding  this  demonstration  we  note  that,  the  polychromatic 
correlation  detection  is  suitable  for  color  coded  signal  detection, 
which  would  be  usable  for  color  coded  pseudorandom  code  application. 

VII.  Conclusion 

This  study  has  led  to  certain  definite  conclusions.  In 
particular: 

1.  Incoherent  optical  correlator  offers  a  higher  signal- 
to-noise  ratio  as  compared  with  a  coherent  correlator. 

2.  The  incoherent  correlator  is  very  suitable  for  decoding 
the  pseudorandom  PN  code  as  applied  in  spread  spectrum 
communication. 

3.  Optimum  matched  filtering  is  a  useful  decoding  tool 
for  burst-error-correlating  PN  code. 

4.  The  technique  is  very  suitable  for  color  signal  correlation 
detection,  which  is  one  step  closer  to  the  actual  human 
recognition, 

5.  The  polychromatic  correlation  technique  is  suitable  for 
color  coded  pseudorandom  code  decoding. 

6.  It  is  feasible  of  utilzing  this  incoherent  correlator 
in  a  real-time  mode,  a  research  program  is  currently 
pursuing. 

7.  The  incoherent  optical  correlator  has  a  superb  noise 
suppression  effect  as  compared  with  the  coherent  system, 
particular  with  the  utilization  of  a  broad  spectral 
band  light  source. 
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